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Preface 



The authors prepared the following booklet in order to make several useful 
topics from the theory of special functions, in particular the spherical har- 
monics and Legendre polynomials of W, available to undergraduates studying 
physics or mathematics. With this audience in mind, nearly all details of the 
calculations and proofs are written out, and extensive background material is 
covered before beginning the main subject matter. The reader is assumed to 
have knowledge of multivariable calculus and linear algebra (especially inner 
product spaces) as well as some level of comfort with reading proofs. 

Literature in this area is scant, and for the undergraduate it is virtually 
nonexistent. To find the development of the spherical harmonics that arise in 
M 3 , physics students can look in almost any text on mathematical methods, 
electrodynamics, or quantum mechanics (see [1], [2], [6], [8], [11], for example), 
and math students can search any book on boundary value problems, PDEs, 
or special functions (see [3], [13], for example). However, the undergraduate 
will have a very difficult time finding accessible material on the corresponding 
topics in arbitrary W. 

The authors used Hochstadt's The Functions of Mathematical Physics [5] 
as a primary reference (which is, unfortunately, out of print). If the reader 
seeks a much more concise treatment of spherical harmonics in an arbitrary 
number of dimensions written at a higher level, [5] is recommended. Much of 
the theory developed below can be found there. 
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Chapter 1 

Introduction and Motivation 



Many important equations in physics involve the Laplace operator, which is 
given by 

d 2 d 2 

2 dx 2 ~^ dy 2 ' ^ ^ 
d 2 d 2 d 2 

3 dx 2 ^ dy 2 ^ dz 2 ' ^ 

in two and three dimensions 1 , respectively. We will see later (Proposition 2.1) 
that the Laplace operator is invariant under a rotation of the coordinate sys- 
tem. Thus, it arises in many physical situations in which there exists spherical 
symmetry, i.e., where physical quantities depend only on the radial distance 
r from some center of symmetry O. For example, the electric potential V in 
free space is found by solving the Laplace equation, 

A$ = 0, (1.3) 

which is rotationally invariant. Also, in quantum mechanics, the wave function 
ij) of a particle in a central field can be found by solving the time-independent 
Schrddinger equation, 

h 2 



2 m A + V ^ ] 



i/; = Etp (1.4) 



where h is Planck's constant, m is the mass of the particle, V(r) is its potential 
energy, and E is its total energy. 

We will give a brief introduction to these problems in two and three di- 
mensions to motivate the main subject of this discussion. In doing this, we 
will get a preview of some of the properties of spherical harmonics — which, 
for now, we can just think of as some special set of functions — that we will 
develop later in the general setting of W. 



x We may drop the subscript if we want to keep the number of dimensions arbitrary. 
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1. Introduction and Motivation 



1.1 Separation of Variables 
Two-Dimensional Case 

Since we are interested in problems with spherical symmetry, let us rewrite the 
Laplace operator in spherical coordinates, which in M 2 are just the ordinary 
polar coordinates 2 , 



-l (V 

i =yj-T y~, (f> — ian 

Alternatively. 



y^+y 2 , = taring). (1.5) 



x = r cos <p, y = r sin <p. 
Using the chain rule, we can rewrite the Laplace operator as 

A 2 ~ — + 1 — + — — (16) 
dr 2 r dr r 2 d(j) 2 

In checking this result, perhaps it is easiest to begin with (1.6) and recover 
(1.1). First we compute 

d dx d dy d , <9 . , d 

= + = COS( ^ +sin0— , (1.7) 

or Or ox or oy ox dy 

which implies that 

d 2 d 2 d 2 d 2 
— -k = cos 2 (p—; + 2 sin (f> cos 4>^-^- + sin 2 (j)-^, 

or* ox* oxoy oy* 

and 

d _ dx d ^ dy d _ . ^ d ^ d 

dip dip dx dip dy dx dy ' 

which gives 

d 2 d 2 d 2 d 2 

u 2 ■ 2 ± u r, 2 ■ i i ° ,2 2 i u 

——, j = r sin 0— r - 2r sin (p cos <p - - + r cos <Ptt^- 
a0 z oar dxdy dy* 

Inserting these into (1.6) gives us back (1.1). 
Thus, (1.3) becomes 

d 2 <5> 1 d§ 1 a 2 $ _ Q 
dr 2 r dr r 2 dip 2 



2 The polar angle cj> actually requires a more elaborate definition, since tan 1 only produces 
angles in the first and fourth quadrants. However, this detail will not concern us here. 
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To solve this equation it is standard to assume that <3?(r, cj>) = x( r )Y((fi), where 
x(f) is a function of r alone and Y{<p) is a function of (j) alone. Then 

y d*x + Y dx + X d*Y = Q 
dr 2 r dr r 2 d(j> 2 

Multiplying by r 2 j\Y and rearranging, 

r 2 d 2 x r dx —1 d 2 Y 
X dr 2 X dr Y d(p 2 

We see that a function of r alone (the left side) is equal to a function of (f> 
alone (the right side). Since we can vary r without changing eft, i.e., without 
changing the right side of the above equation, it must be that the left side of 
the above equation does not vary with r either. This means the left side of the 
above equation is not really a function of r but a constant. As a consequence, 
the right side is the same constant. Thus, for some —A we can write 

— 1 d?Y ^ r 2 d 2 x dx ^ ^ 

Y d(j) 2 x dr 2 x dr 

We solve Y" = XY to get the linearly independent solutions 



Y{4>) 



e^, e _v ^ if A > 0, 

1,0 if A = 0, 

sin(v^), cos(vlAfa) if A < 0, 



but we must reject some of these solutions. Since (ro, 4>o) represents the same 
point as (ro, (fio + 2irk) for any k G Z, we require Y(cfi) to have period 2ir. Thus, 
we can only accept the linearly independent periodic solutions 1, sin ^\/jA|V>^ , 

and cos f\/|A[0j , where \/|A| must be an integer. Then, let us replace A with 
— m 2 and write our linearly independent solutions to Y" = —m 2 Y as 

Y\ t n = cos (n<f>), Y 2y m = sin (m4>), (1.9) 

where 3 n G No and m G N. Notice from (1.6) that d 2 /d(j) 2 is the angular 
part of the Laplace operator in two dimensions and that the solutions given in 
(1.9) are eigenfunctions of the d 2 /d(j) 2 operator. We will see in Chapter 4 that 
the functions in (1.9) are actually spherical harmonics; however, since we have 
not yet given a definition of a spherical harmonic, for now we will just refer to 
these as functions Y . The reader should keep in mind that characteristics of 
the Y's we comment on here will generalize when we move to W . 



3 We use the notation N = {1, 2, . . . } and No = {0, 1,2,...}. 



4 



1. Introduction and Motivation 



We can also solve easily for the functions x( r ) that satisfy (1.8), but this 
does not concern us here. Let us instead notice a few properties of the functions 
Y. Let us consider these to be functions r m sin (mcfi), r n cos (n</>) on R 2 that 
have been restricted to the unit circle, where r = 1, and let us analyze the 
extended functions on R 2 . We will first rewrite them using Euler's formula, 
e 1 ^ = cos cp + i sin cf>, which implies 

(x + iy) n = (re^Y = r n e in(p = r n [cos (ncj)) + isin (ncj))\ , 
(x - iy) n = (re-^Y = r^e"^ = r n [cos (wf>) -ism (n<£)] , 

so that 

r n cos (n4>) = -[{x + iy) n + (x - iy) n ] = f H ljti (x, y), 
r n sin (ru/>) = ^ [(z + iy) n - (x - iy) n ] = f H 2 , n {x, y). 

We notice that the Y's can be written as polynomials restricted to the unit 
circle, where r = 1. Furthermore, observe that 

Hi >n (tx, ty) = t n H ltU (x, y), and H 2 , n (tx, ty) = t n H 2>n {x, y); 

we call polynomials with this property homogeneous of degree n. Moreover, 
the reader can also check that these polynomials satisfy the Laplace equation 
(1.3), by either using (1.1) or the Laplace operator in polar coordinates (1.6) 
for the computation, i.e., 

A 2 fli,„ = 0, and A 2 H 2 , n = 0. 

Let us also notice that the y's of different degree are orthogonal over the unit 
circle, which means 

Jq 77 sin (ncf)) sin (m(f>) d<p = 0, if n ^ m, 
Jq W cos (ncf)) cos (m<j)) dcfi = 0, if n ^ m, 
Jq 77 sin (ncf)) cos (mcfi) dcj> = 0, if n ^ m, 

as we can easily compute by taking advantage of Euler's formula. For instance, 
we can calculate 
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or 



which becomes 

i(m+n)<f> g— i(m+n)<j> gi(m—n)<l> _|_ ^—i{m—n)<j> 

2 ~~ 2 

y ^ sin[(m + n)<^>] — cos[(m — n)0]^ d</> = 0, 
o 

for n ^ m. The reader can check the rest in a similar fashion. 

Finally, by recalling the theorems of Fourier analysis 4 , we know that any 
"reasonable" function defined on the unit circle can be expanded in a Fourier 
series. That is, given a function / : [0, 2ir) — > M satisfying certain conditions 
(that do not concern us in this introduction) , we can write 

oo oo 

f{4>) = ^ a,m sin ( m ^) + bn cos ( n( ^ 

m=l n=0 

for some constants a m ,b n . We say that the Y's make up a complete set of 
functions over the unit circle, since we can expand any nice function f(4>) 
defined on [0, 2ir) in terms of them. 

Before we move on, we will show that (1.4) can be approached using a 
method almost identical to the method we used above. Let us assume that 
the solution -0 to (1.4) can be written as ijj(r,(f>) = x(r)Y(<^>). Then, using 
polar coordinates, the equation becomes 

Y (- — f- -Jl—^ + V (r)] x+— (- — —) Y - E X Y. 
\ 2m dr 2 2mr dr J r 2 \ 2m d(j) 2 J 

Multiplying by r 2 /xY and rearranging, 

X V 2m dr 2 2m dr J Y \ 2m dip 2 J 

Once again, we see that a function of r alone is equal to a function of (ft alone, 
and we conclude that both sides of the above equation must be equal to the 
same constant. Using the same reasoning as before, we write this constant 
as — £ 2 h 2 /2m, where £ is an integer. If we carry out the calculation, we will 
see that the functions Y(<p) are the same ones we found previously in this 
subsection. We will also find the radial equation 

1 / r 2 h 2 d 2 rh 2 d ~ T , A ^ 2 i 2 n 2 
-[- — —, -— 1Z+ r 2 V(r) )x-Er 2 - 



2m dr 2 2m dr 1 2m ' 



4 Doing so will not be necessary to understand the material we present here, but the reader 
unfamiliar with Fourier analysis may choose to consult [3]. 
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1. Introduction and Motivation 



which we can rewrite as 



2m \ dr 2 



ld_ 

r dr 



+ 



V(r) + 



l h 2 



2mr z 



X = Ex- 



It is interesting to note that this equation resembles (1.4). If we think of r 
as our single independent variable and x as our wave function, we have an 
effective potential energy 



V eS (r) = V(r) + 



l h 2 



2mr 2 



where we call the second term the centrifugal term. Thus x represents a 
fictitious particle that feels an effective forece 



cff 



cff 



dr 



r. 



We see that the centrifugal term contributes a force 



F. 



centrifugal 



£ 2 h 2 

3 r ' 

mr* 



pushing the particle away from the center of symmetry. 



Three-Dimensional Case 

Here, we will follow a procedure almost identical to that of the last subsection, 
but we will not take the discussion as far in three dimensions. The Y's we will 
find in three dimensions are more widely used than those in any other number 
of dimensions; however, a thorough development of the functions in M 3 could 
take many pages, and this would distract us from our goal of moving to p 
dimensions. Furthermore, the results that we would find in three dimensions 
are only special cases of more general theorems we will develop later in the 
discussion. If the reader is interested in studying the usual spherical harmonics 
of M 3 in depth, there are a multitude of sources we can recommend, including 
[1], [2], [3] and [13]. 

Let us rewrite the Laplace operator in spherical coordinates, which are 
given by 



r = y / x 2 + y 2 + z 2 , 6 = tan" 1 + y j ; = tan" 1 (|) . (1.10) 

Alternatively, 



x = r sin 8 cos eft, y = r sin 8 sin 0, z = rcosd. 
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Using the chain rule, we find 



A, = ~ 



1 8 



Or 



Or 



+ 1 



1 8 
sin 6 86 



sin 



86 



+ 



8 2 



sin 



(1.11) 



As in the two-dimensional case, it is probably easiest to verify this formula by 
starting with (1.11) and producing (1.2). The reader should check this result 
for practice, proceeding exactly as we did beginning in (1.7). 
Inserting (1.11) into (1.3) gives 



r 2 8r 



,0$ 

Or 



+ 



1 



1 8 
sin 6 86 



sin ( 



} 0$ 
50 



+ 



1 <9 2 $ 



sin 



0. 



Searching for solutions <3? in the form <3?(r, 0, 0) = x( r )Y(6, (f>), where x( r ) is a 
function of r alone and Y(6, 4>) is a function of only and <f>, this becomes 



y 



1 _d_ 

dr 



dr 



* + ~2 



1 



sinfldfl V Sm 50 



9 



1 9 



2 1 



sin 



Y = 0. 



Multiplying by r 2 /yY and rearranging, 



1 _d_ 



dr 



X 



-1 



i_ a_ 



sin 



sin 



+ 



i a : 



.2 1 



sin 



y. 



Since we have found that a function of r alone is equal to a function of only 
6 and </>, we use the same reasoning as in the previous subsection to conclude 
that both sides of the above equation must be equal to the same constant, call 
it —A. This implies that 



1 8 

sin 6 86 



sin 



,1 

86 



+ 



1 8 2 



sin 



Y = XY. 



(1.12) 



At this point, we will stop. It turns out that the Y's which satisfy this equation 
are actually spherical harmonics. Comparing (1.11) and (1.12), we see that 
the y's are eigenfunctions of the angular part of the Laplace operator, just 
as in two dimensions. If we studied the functions Y further we would also 
find that, analogously to what we noticed in 1R 2 , the y's form a complete set 
over the unit sphere, each Y is a homogeneous polynomial restricted to the 
unit sphere, and these polynomials satisfy the Laplace equation. However, 
to develop these results and the many others that exist would require us to 
study Legendre's equation, Legendre polynomials, and associated Legendre 
functions, and we will choose to leave such an in-depth analysis to the general 
case of p dimensions. 

We will now move on to see how these functions Y are related to angular 
momentum in quantum mechanics. 
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1. Introduction and Motivation 



1.2 Quantum Mechanical Angular Momentum 

We have seen that spherical harmonics in two and three dimensions relate to 
the one-dimensional sphere (circle) and the two-dimensional sphere (surface 
of a regular ball) respectively. In quantum mechanics, rotations of a system 
are generated by the angular momentum operator. Spherical symmetry means 
invariance under all such rotations. Therefore, a relation between the theory of 
spherical harmonics and the theory of angular momentum is not only expected 
but is a natural and fundamental result. 

Recall that in classical mechanics, the angular momentum of a particle is 
defined by the cross product 

L = r x p, 

where p is its linear momentum. To find the quantum mechanical angular 
momentum operator, we make the substitution pi i— >■ —ihd/dxi where x\ = x, 
X2 = y, and X3 = z. Thus, we see that 5 



L = —ihfx V 

where 

V : 



Odd 
dx ' dy ' dz 



Two-Dimensional Case 



In the plane, the angular momentum operator has only one component, given 
by 

L = -ih x- y— 

\ ay ox 

Using the polar coordinates defined in (1.5) and the chain rule, we can rewrite 
this as 

i = ~ ih h 

as the reader should have no problem checking using the same strategy s/he 
used to verify (1.6) and (1.11). Then 

L ~ n dp* 

and we check using (1.9) that the functions Y(<f>) are eigenfunctions of the L 2 
operator. In particular, 

L 2 Y mtj {4>) = -h 2 ^Y mJ = h 2 m 2 Y mJ , 



5 The hat above the angular momentum indicates that it is an operator (not a unit vector) 
in this section. 
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and we see that the function Y m j is associated with the eigenvalue h 2 m 2 . 

In quantum mechanics, operators such as L represent dynamical variables. 
If an operator O has eigenfunctions ipk with corresponding eigenvalues 
then a particle in state ipk wu l be observed to have a value of for the 
dynamical variable O. Therefore, we see that a particle in the state Y m j will 
be observed to have a value of h 2 m 2 for its angular momentum squared. We 
say the function Y m j carries angular momentum Km. 



Three-Dimensional Case 

Things work similarly in three dimensions, where the angular momentum op- 
erator has the components 

L x = -ih y— - z— , 
\ oz ay ) 

L z = -%h x- y— 

\ ay ox 

Using the spherical coordinates defined in (1.10) and the chain rule, we 
can rewrite these as 

/ d d 
L x = ih\ sin 9— + cos d> cot 9 — 
\ 06 da 



t ( n 9 . , n d 

L v = —in cos fl— + sin d> cot 9 — 
\ 69 oq 

o<p 



the reader should check these formulas. These equations allow us to compute 

-2 



L 2 = L ■ L = L 2 X + L 2 + Li. Carrying out the multiplication, 



L 2 = -h 2 



l d ( n d\ l & 

smo— + 



sin 9 89 V 89 sin 2 



)2 



(1.13) 



and we see that by (1.12), the functions Y are eigenfunctions of the L 2 oper- 
ator. 

We claimed in the last section that the functions Y{9, 4>) were homogeneous 
polynomials with restricted domain, so let us write Yi(9,(p) where £ denotes 
the degree of homogeneity. In Section 4.2, we will see an easy way to compute 

the eigenvalue of L 2 associated with Ye, and it will turn out to be h 2 £(£ + 1). 
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1. Introduction and Motivation 



So we claim that in three dimensions, the function Yi(9, 4>) carries an angular 
momentum of %\J l{l + 1). 

In Chapter 4, we will give rigorous foundations to the seemingly coinci- 
dental facts we have discovered in this chapter about the functions Y that 
arose as solutions to certain differential equations. But first, we will devote a 
chapter to gaining some practice and intuition working in M p . 



Chapter 2 

Working in p Dimensions 



In this chapter, we spend some time developing our skills in performing calcu- 
lations in MP and exercising our abilities in visualizing a p-dimensional space 
for an arbitrary natural number p. We will use the majority of the results 
we obtain here in the development of our main subject, but some topics we 
discuss just out of pure interest or to improve our intuition. 

First, let us generalize the definition of the Laplace operator to M p , where 
a point 1 x is given by the ordered pair (x\, X2, ■ ■ ■ , x p ). 

Definition The Laplace operator in MP is given by 



Let us quickly consider orthogonal rotations of the coordinate axes in MP. 
Such rotations leave the lengths of vectors unchanged. Indeed, the length of 
a vector is a geometric quantity; rotating the coordinate system we use to 
describe the vector leaves its length invariant. In fact, in a more abstract 
setting, we could define a rotation to be any transformation of coordinates 
that leaves the lengths of vectors unchanged. 

In what follows, we let x denote a column vector 2 (x\,X2, ■ ■ ■ ,x P Y in MP 
and use (•,•) to represent the dot product of two vectors. The fact that a 

x We will not place vector arrows above points x in R p . 

2 The superscript t denotes the operation of matrix transposition. 




(2.1) 




2.1 Rotations in W 



n 
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2. Working in p Dimensions 



rotation matrix R leaves the length of x invariant means {Rx,Rx) = (x,x). 
Moreover, since the dot product between any two vectors x, y can be written 
as 

v) = 2 i( x + y> x + y) ~ ( x > x ) ~ (y> 2/)) > 

it follows that coordinate rotations leave all dot products invariant. 

Notice further that we can write dot products such as (x, y) as matrix prod- 
ucts y t x. In this notation, the requirement that (Rx, Ry) = (x, y) translates 
into the necessity of {Ry) t {Rx) = y t x, or y t R t Rx = y l x. Since this equation 
must hold for all x, y G we can conclude that any rotation matrix R must 
satisfy R l R = I, where / is the identity matrix. 

Now we can verify our claim in the first few sentences of this booklet that 
the Laplace operator A p remains unchanged after being subjected to a rotation 
of coordinates. 

Proposition 2.1 The Laplace operator A p is invariant under coordinate ro- 
tations. That is, if R is a rotation matrix and x' = Rx, then A' = A 



■p, 



i.e. 




E 



d 

dxj 



Proof This can be proved very easily by noticing that A p = V p • V p is a 
dot product of vector operators. Since all dot products are unchanged by 
coordinate rotations, we can conclude that A p is not affected by any rotation 
R. 

In case the reader is not satisfied with this quick justification, let us com- 
pute Ap, the Laplace operator after application of a rotation of coordinates 
R. Since R is a rotation matrix, it is orthogonal, i.e. RR l = I. Then, using 
the chain rule, 



A, 



E 

3=1 
V 

£ 

3=1 



d 

dxj 



^2 Rkj dx> 

\k=i k 




ydx^A_ 

dx 3 dx 'k 



V) 



d 
dx'n 



^ dx' e d 
d Xj dd t 



so 



A 



_ >p d d 




p 

E 

fc=i 



d 
dx',„ 



A'. 



This proves the proposition. 
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2.2 Spherical Coordinates in p Dimensions 

Now, in order to develop some experience and intuition working in higher- 
dimensional spaces, we will develop the spherical coordinate system for W in 
considerable detail. In particular, we will use an inductive technique to come 
up with the expression of the spherical coordinates in p dimensions in terms 
of the corresponding Cartesian coordinates. 

We will let our space have axes denoted x\, x 2 , ■ ■ ■ First, in two dimensions, 
spherical coordinates are just the polar coordinates given in (1.5), 

r = \Jx\ + x\ G [0, do), 
(ft = tan _1 (x 2 /xi) G [0,271"), 

where r is the distance from the origin and (ft is the azimuthal angle 3 in the 
plane that measures the rotation around the origin. The inverse transforma- 
tion is 

x\ = r cos 4> , X2 = r sin (ft . 

When we move to three dimensions we add an axis, naming it x$, per- 
pendicular to the plane. Now, the polar coordinates above can only define a 
location in the plane; thus, they only tell us on which vertical line (i.e., line 
parallel to the X3-axis) we lie, as we can see in Figure 2.1 with p = 3. To 
pinpoint our location on this line, we introduce a new angle Q\. When we also 
redefine r to be the three-dimensional distance from the origin, we have the 
spherical coordinates given in (1.10), 

r = v / xf~+xf+lcf G [0, oo), 
(ft = tan" 1 (x 2 /xi) G [0,271"), 

9i = tan -1 (yjx\ + x\/x^ G [0,7r]. 

Now let us imagine moving to four dimensions by adding an axis — name 
it the X4-axis — perpendicular to the three-dimensional space just discussed. 
The 3D spherical coordinates given above can only define a location in M 3 , so 
they only tell us on which "vertical" line (i.e., line parallel to the X4-axis) we lie, 
as in Figure 2.1 with p = 4. We thus introduce a new angle Q% to determine the 
location on this line. We redefine r to be the four-dimensional distance from 

3 We keep the definition of 4> sloppy throughout this section. A more precise formula 
would use a two-argument tan -1 function that produces angles on the entire unit circle. 
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2. Working in p Dimensions 




Figure 2.1: In going from W^ 1 to W, we visualize R p_1 as a plane and add a new 
perpendicular direction. We introduce a new angular coordinate 8 p -2 to determine the 
location in the new direction. 



the origin, and this completes the construction of the 4D spherical coordinates, 



0-2 



yjx\ + x\ + X% + x\ 

tan -1 (a^/zi) 
tan -1 



G [0,oo), 
G [0,2vr), 

G [0,vr], 



! x\ + xl/x 3 
tan -1 ^ y 7 ^! + x l + x i/ x 4^ G [0,tt], 



where it should be clear from the figure that the new angle 62 only ranges 
from to 7T. 

In going from three to four dimensions, we mimicked the way in which we 
transitioned from two to three dimensions. We follow the same procedure each 
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time we move up in dimension. For arbitrary p, this yields 



rv>2 i ™2 i ... _i ™2 



tan 1 (X2/X1), 



tan 1 ( + x \l ^3 



6^-2 = tan 1 Wzf + ;r| H ^ x 1-\j r - 

where the ranges on the coordinates are as expected from the previous cases. 

We have thus defined the spherical coordinates in p-dimensions in terms of 
the corresponding Cartesian coordinates. To write down the inverse relations 
we use projections, with Figure 2.1 as an aid. As before, we will derive these 
relations in detail for a few instructive cases before writing down the most 
general expressions. 

We have already written down x\, X2 in terms of r, (ft so let's use M 3 as our 
first example. We imagine M 3 as the direct sum of a two-dimensional plane M 2 
with the real line M. Then, given the point P in IR 3 with spherical coordinates 
(r3, (ft, 9i), the vector ofi can be written as a sum OA + 0$, where OA lies 
along the X3-axis and has magnitude r% cos#i and OB^ lies in the plane and 
has magnitude T2 = r% sin#i. The point B thus has spherical coordinates 
(r2,(p) in the plane, implying 

xi = T2 cos 4> and X2 = r2 sin (ft, 

so 

x\ = r% sin 0i cos (ft , (2-2) 

X2 = ^3 sin#i sin^ , (2-3) 

£3 = ^3 cos Q\ . (2-4) 

Let us examine the situation in M 4 using the same technique. Given the 
point P with radial distance from the origin, we decompose 0~P into two 
vectors OA+OB^, where OA lies along the X4-axis and has magnitude cos O2 
and OB lies in the "plane" and has magnitude r% = sin 62- The point B 
thus has spherical coordinates (rs,(ft, #i) in the three-dimensional space, so 
x%,X2,X3 are as in (2.2)-(2.4). Therefore 

x\ = r4 sin 82 sin 9\ cos (ft , 
X2 = r A sin #2 sin#i sm.<ft , 
X3 = r4 sin #2 cos 9\ , 

X4 = T4 COS 6*2 • 
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2. Working in p Dimensions 



If the reader has understood the previous constructions, the expressions 
in p dimensions should be evident. Given the radius r p in IR P , we project the 
position vector onto M p_1 , obtaining the radius in the subspace given by 
r p _i = r p sin# p _2. In this way we can perform a series of projections to get 
down to a space in which we already know the relations. This procedure leads 
to (dropping the subscript on r p ) 

x\ = r sin# p _2 sin#p_3 • • - sin 63 sin #2 sin#i cos</>, 

X2 = r sin 9 p -2 sin 9 p -3 ■ ■ • sin 63 sin #2 sin 9\ sin <j>, 

X3 = r sin6'p_2 sin 6^-3 • • • sin #3 sin 02 cos#i, 

X4 = r sin #p_2 sin 6 p -3 ■ ■ ■ sin 83 cos #2, 

Xp-i = r sin 6* p _ 2 cos Q p - 3, 
x p = r cos9 p -2- 

With this chapter as an exception, we will rarely refer explicitly to the 
angles (f>, #1, . . . , Q p -2- However, we will frequently use r = \J x\ + • • • + x|. 

2.3 The Sphere in Higher Dimensions 

We will give definitions of the sphere and the ball in an arbitrary number of 
dimensions that are analogous to the definitions of the familiar sphere and ball 
we visualize embedded in M 3 . 

Definition The (p — l)-sphere of radius 5 centered at xq is the set 

S^ixo) = {x G R p : \x - xq\ = 5}. 

The unit (p — l)-sphere A is the set S p ~ l = f S p l (0). 

Definition The open p-ball of radius 5 centered at x$ is the set 

B P (x ) = {xeW :\x-x \ <5}. 

The open unit p-ball is the set B p = f B p (0). The closed p-ball is the set 
BP = B p US p - 1 . 

Notice we call the sphere that we picture embedded in W the (p — 1)- 
sphere because it is (p — l)-dimensional. It requires p — 1 angles to define 
one's location on the sphere, as we saw in Section 2.2. However, it requires p 

4 Frequently, we will drop the "unit," though it is still implied. 
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coordinates to locate a point in the ball because we must also specify r; this 
justifies the notation B p . Notice also that the (p — l)-sphere is the boundary 
of the p-ball when we think of these as subsets of W; we write = dB p . 

Let us compute the surface area of the (p — l)-sphere. Towards this goal, 
we recall that that the gamma function is defined by 



r(z) = J e-H"- 1 dt, 



for any z € C such that Re(z) > 0. Then 

Lemma 2.2 For all p G C such that Re(p) > 0, we have 



rOO 



/ 



1 /p 



e ■ rP- 1 dr = - F . 

2 V2 



o 

2 



Proof Using the substitution u = r 

+oo oo +oo 

e-V" 1 dr = [ = 1 I e^u^ du= l v( P 

J 2^u 2 J 2 V2 



as sought. ■ 

Proposition 2.3 // denotes the solid angle in MP (equivalent numeri- 
cally to the surface area of S p ~ l ), then 

ilp - l ~T{p/2y 



Before proving this we will digress slightly and discuss how the surface area 
of the (p — l)-sphere relates to the volume of the p-ball. First, consider the 
p-ball of radius r centered at the origin. Notice that the radius r completely 
determines such a ball. If we were to determine the volume V p of this p-ball, we 
would obtain V p = cr p where c is some constant. Here, we have determined 
that V p oc r p because V p must have dimensions of [length] p and the only 
variable characterizing the p-ball is r (which has dimensions of length). Now 
if we differentiate V p with respect to r, we get 

d Yp = (p-i) cr p -\ 
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2. Working in p Dimensions 



Thus, if the radius of a p-ball changes by an infinitesimal amount dr, its volume 
will change by some infinitesimal amount 6V P , and 

SV P = (p- l)cr p ~ l Sr. (2.5) 

But in this case, the small change in volume 5V should equal the surface area 
A p -i(r) of the p-ball multiplied by the small change in radius 5r, 

6V p = A p _ 1 (r)Sr. (2.6) 

From (2.5) and (2.6) we see that 

A p „ l {r) = {p-l)crV-\ 

and if we let f2 p _i = (p — 1) c denote the numerical value of the surface area 
when r = 1 (as in the statement of Proposition 2.3), we get 5 

A p _i(r) =n^ ir P~ 1 . (2.7) 

We see that if we want to carry out an integral over W when the integrand 
depends only on r, we can use the differential volume element 

dV p = Ap-i(r) dr = r p_1 Op_i dr. 

Now we are ready to prove Proposition 2.3. 

Proof Consider the integral, 

oo oo oo 

J= j dx! J dx 2 - - J dx p e-^ +x " + - +x p\ 

— oo — oo — oo 

This is really 

/CO \ P 

J= J J e~ x2 dx J = (y^) p . 
Using spherical coordinates however, we can write 



J= J J dV p e~ r2 = fip_i J e" r ' 2 r^ 1 dr, 
sp- 1 



5 We should stress that, since a physicist assigns physical dimensions to each and every 
quantity, s/he would differentiate between the solid angle fi p _i and the surface area A p -i of 
the (p— l)-sphere. Although numerically the two quantities are equal for the unit sphere since 
r — 1, they are different quantities since f2 p _i is dimensionless while A p -i has dimensions 
of [length] P" 1 . 
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since O p _i is just a constant. Therefore 



n 



p-i 



jco ^ r 2 . r p 1 

and, with the help of Lemma 2.2 we arrive at the advertised result. 



Remark As a check, we can use this formula to determine the surface area 
(or circumference) of the 1-sphere (or circle) as well as the surface area of the 
2-sphere, which is the familiar sphere that we embed in M 3 . As expected, 

2vr 3 / 2 

^ = 2*. ^ = ^ = 4*, 

using T(3/2) = y/ii/2. It is interesting to consider the 0-sphere, i.e., the sphere 
that we visualize embedded in M.. S° consists of all points on the real line that 
are unit distance from the origin, so 5° = {—1, 1}. Using Lemma 2.3, we find 
£Iq = 2, since T(l/2) = y/n. That is, the surface area of just two points on the 
real line is finite! This is consistent with standard calculus, though. On the 
real line, the radial distance is the absolute value of x. Hence, the concept of 
a function f(r) possessing spherical symmetry coincides with the concept of 
an even function, f(x) = f(—x). Then 

+oo oo oo 

J f(x)dx = 2 J f(x)dx = J f(r)n r°dr. 

-oo 



2.4 Arc Length in Spherical Coordinates 

Now, let us compute the formula for the infinitesimal arc length 5s associated 
with a small displacement along a curve in M p . In this section, we use infinites- 
imals in the physicist's way; that is, 5 s and the other small lengths involved 
represent tiny displacements that are small enough so that the errors involved 
in our approximations (such as assuming a small portion of a curve is almost 
straight) are as small as desired (say, less than some given e > 0). We only use 
the formulas we derive under this assumption in the limit as 5s — > and like- 
wise with the other infinitesimal quantities, so our approximations introduce 
no error into future calculations. 

Working in M p , let us use the spherical coordinates we developed in Section 
2.2. Consider a point P with coordinates (r, (f>, 9%, . . . , 6 p -2)- Consider a small 

displacement vector 5s = PP' at the point P in an arbitrary direction. One 
component of this displacement is along the radial direction and has magnitude 
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2. Working in p Dimensions 



5r. The other p — 1 components of ets are orthogonal to 5f and thus 6 are 
along the surface 5r _1 (0). Thus by the Pythagorean Theorem, we have so far 
5s 2 = 5r 2 + • • • where the rest to be filled in involves small changes in the 
angles 6» p _ 2 , . . . ,0i,(f>. 




O 




5r P' 



Figure 2.2: The plane containing the points 0,P,P'. The displacement vector Ss has a 
radial component of magnitude Sr and a component along the (p — l)-sphere of magnitude 
r 8ojp-i. 



Let us consider the plane containing the vectors OP and PP', where O is 
the origin (see Figure 2.2). We are familiar with the formula for an infinitesimal 
arc length in this two-dimensional plane, namely 

Ss 2 = 5r 2 + r 2 5u 2 l: (2.8) 



where 5u p -\ is the small angle between the vectors OP and OP'. Note 5uo p -\ 
is the angular displacement that occurs along the unit (p — l)-sphere. 

Let us constrain the displacement 5s to lie along a sphere by requiring 
5r = 0. If we compute 5s in this special case, we have found 5u} p -i = 5s /r 
which we can plug into (2.8) to obtain the general result. 

We can determine 5u p -i inductively, starting in the plane with p = 2. We 

6 The remaining components of 5s lie in a surface where Sr = or r = const. This is a 
[p — l)-sphere. 
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know that on the unit circle 5s = r 5(f), so 

5^ = ^ = 5cf. 

In M 3 , 5s has a component r 50\ along the direction of increasing 0%, so that 
5s 2 = r 2 50 2 or <5u;| = 50 2 + • • • , and the remaining orthogonal components are 
parallel to the M 2 plane. Projecting PP' onto M 2 , we see that the infinitesimal 
displacement lies on a circle of radius sin 0\ . Using our result from M 2 , we get 

5u\ = 50\ + sin 2 0i 5oo 2 = 56\ + sin 2 1 5(j) 2 , 

the familiar result in M 3 . 

We can move to M in an analogous fashion. An infinitesimal arc length 
5s on the 3-sphere has a component r 502 along the direction of increasing 
02, so that 5s 2 = r 2 50 2 or 5u 2 = 50\ + • • • , and the remaining orthogonal 

components are parallel to M 3 which we picture as a plane. Projecting PP' 
onto M 3 , we see that the infinitesimal displacement lies on a 2-sphere of radius 
sin #2 , so 

5u% = 561 + sin 2 9 2 5ujI = 50\ + sin 2 P ^ 2 S0j + sin 2 2 sin 2 1 5(f) 2 . 

It is clear that this pattern continues, so that in R p , 

5u} p -i = 50 2 _ 2 + sin 2 p -2 5uj 2 _ 2 . 

By filling in the form of 5lo 2 _ 2 in this expression, we get 

5uj p -i = 50 2 _ 2 + sin 2 6»p„ 2 (50 2 _ 3 + sin 2 P ^ 3 (5w 2 _ 3 ). 

Continuing to expand this expression, we come up with 

5uJ p -i = 50 p _ 2 + sin p -2 $0 p -3 + sin 6 p -2 sin p - 3 50 p _ 4 
+ • • • + sin 2 p -2 sin 2 p s ■ ■ ■ sin 2 0\ 5cj) 2 . 

Putting this into (2.8) completes our computation. 

2.5 The Divergence Theorem in W 

Let us now return to the use of Cartesian coordinates to give an intuitive 
"derivation" of the divergence theorem in p dimensions, which we will use 
numerous times in Chapter 4. For brevity, our justification of this theorem 
will be rather physical. For a more mathematical treatment, the interested 
reader should consult a calculus book (for example, [12]) for the theorem in 
M 3 and see if s/he can generalize the proof there to W. For a rigorous proof 
of the divergence theorem in p dimensions, s/he may consult an analysis text 
(for example, [10]). 
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2. Working in p Dimensions 



Theorem 2.4 Let F be a continuously differentiable vector field defined in 
the neighborhood of some closed, bounded domain V in W which has smooth 
boundary dV . Then 1 

I V ■ FdPx = j F -nda, (2.9) 

V dV 

where h is the unit outward normal vector on dV and da is the differential 
element of surface area on dV . 

"Proof" We will interpret F as the flux density of some p-dimensional fluid 
moving through the volume V. In unit time at a point x, the volume of fluid 
which flows past an arbitrarily oriented unit surface with unit normal vector 
h is given by F(x) ■ n. 

Let us fill the interior of V with a "grid" of disjoint boxes, none intersecting 
the boundary dV — see Figure 2.3. If the boxes are comparable in size to V, 
they will not make a complete covering of the interior, since many regions of 
V near dV will remain uncovered; however, as the lengths of the box edges 
approach zero, the entire interior of V can be covered by these boxes. We 




Figure 2.3: An example of a domain V in M. 2 . A similar picture would describe 
the situation in arbitrary M p . We consider one box Q with bottom-left corner 
at (xi,X2) and top-right corner at (x\ + 8x\,X2 + 8x2)- 

will consider one small box Q in the interior of V, say with one corner at 
(xi,X2, ■ ■ ■ ,x p ) and another at {x\ + 8x\,X2 + 8x2, ...,x p + 5x p ), where the 



Here and from now on, d v x = dxidx2 ■ ■ ■ dx p . 
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all positive. Consider the integral 



J = j> F ■ hda. 

8Q 



The boundary dQ consists of 2p planes, two orthogonal to each coordinate 
axis. The two planes orthogonal to the x\ axis have surface area equal to 
5x2 8x3 • • • Sx p . Since F is continuous, the integral of F ■ h over these surfaces 
can be written as 8 



F(x\ + 8x\,X2-, • • • , x p ) • x\ 8x2 8x3 ■ ■ ■ 5x 



and 



F(x 1 ,x 2 , . . . , x p ) • (-xi) 5x2 8x 3 - ■■ 8x p , 

since F changes very little over this plane as 6x\, 8x2, ■ ■ ■ , 8x p — > 0. Now, 
there was nothing special about the xi-axis, so we see that for each 
the integral J will include a term 

„ 8x1 ■ ■ ■ Sx p 



F(x\ , . . . , X\ -\- 8xi, . . • , Xp) F{x\ , .... x'j, ... , Xp) 

This can be rewritten as 

F[x\ , . . . , xi + 8xi, . . . , Xp) F{x\ , . . . , Xj, . . . , Xp 



which becomes 



dF(xi, . . . ,Xp 
dxi 



8xi 



Xi 8x\ ■ ■ ■ 8x„ or 



8xi 



■ Xi 8x1 ■ ■ ■ 8x p , 



dFj{xi, ...,x p 
dxi 



8X1 ■ ■ ■ 8x r 



in the limit with which we are concerned. Putting all these terms together, 
we get 



dFi(xi, . . . ,Xp) 



i=l 



dxi 



8x\ ■ ■ ■ 8x p . 



But this is just 



V p • F{x\, . . . ,x p ) 8x± ■ ■ ■ 8x p = / V-Fd p x 



Q 



in our limit. We have thus shown that 



V-Fd p x= <b F-hda 

8Q 



3 Here, xi is a unit vector in the direction of increasing x\. 
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2. Working in p Dimensions 



for the special case of an infinitesimal box Q inside V. Let us now sum this 
result over all the boxes inside V, 

Jv-Fd p x = j>F-hda, (2.10) 

Q* q % & dQi 

where we are concerned with the limit as the lengths of the box edges go 
to zero and the number of boxes in the covering approaches infinity. In this 
limit, the volume of the uncovered regions of V near dV approaches zero. 
Since F is continuously differentiable, it follows that V • F is continuous and 
thus bounded over the closed and bounded region V. Thus, the integral of 
V • F over regions in V not covered by boxes approaches zero. The left side 
of (2.10) therefore becomes 



^Jv-FdPx — ► Jv-Fd p x (2.11) 



in this limit. Now, in the right side of (2.10), let us consider all the planes 
which bound the boxes Qi that are included in the sum. Each "interior" plane 
appears twice in the sum, once as the "right" side of one box and a second 
time as the "left" side of another box. In each of these appearances, both F 
and da remain the same but the vector n shows up with opposite sign. Thus, 
the integral over all the interior planes vanishes. The only terms that are not 
canceled in the integral on the right side of (2.10) make up the integral over 
the "exterior" planes, which we write as d{JQi. That is, 



F ■ hda = f F • n da, 

dQr d\JQi 



so using (2.11), 



V-Fd p x= j> F-nda. (2.12) 

V d\JQi 

Now we are very close to (2.9), but there is one difficulty due to the fact 
that the outward normal vector n for d |J i Qi always points along one of the 
coordinate axes while n can point in any direction for dV. To reconcile this 
difference, we realize that integrating F ■ n over a closed surface gives us the 
volume of fluid that has passed through this surface in unit time. The same 
volume of fluid must pass through both d \J { Qi and dV since d |J i Qi — > dV . 
Therefore, 

F ■ hda = j> F ■ h da. 
aUQi dv 
Combining this with (2.12) completes the proof. ■ 
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2.6 A p in Spherical Coordinates 

To compute A p in spherical coordinates, we could use the chain rule. This 
would involve converting all the derivatives with respect to Xj in (2.1) to 
derivatives with respect to spherical coordinates as we did for A2 and A3 in 
Section 1.1. Such an undertaking would be quite messy, however, and adds 
little additional insight. 

Another approach would be to use the general formula from differential 
geometry 

where g^ u is the metric tensor of the space of interest, g = detfc^] and g^ v 
the inverse of the metric tensor. Since we have avoided this route thus far, we 
will not use it here either. 

We can can actually use an integration trick to determine the form of 
the p-dimensional Laplace operator A p in spherical coordinates. This is the 
approach followed below. 



V 



O 




Figure 2.4: A region V in W. 

Consider the p-dimensional cone-like region depicted in Figure 2.4. Let T> 
be the shaded region (a truncated cone), which has left and right boundaries 
given by portions of Sf~(0) and S^7 Sr (Q) respectively. We are concerned with 
the limit where the quantity 5r — > 0. This region spans a solid angle of ui, 
which means the surface area of the region is a fraction uj/£I p -i of the surface 
area of a complete sphere of the same radius. Using (2.7), this implies that 
the surface area of the left boundary is wr p_1 while that of the right boundary 
is u(r + 5r) p ~ l . 

Now let / : W — > 1R be a twice continuously differentiable function, and 
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2. Working in p Dimensions 



use the divergence theorem in p dimensions to find 9 

jA p fdPx = Jv p -(V p f)dPx = Jv p f-nda = J^da, 

V V dV dV 

where h is the external unit normal vector to T> and da is a differential element 
of surface area. By the mean value theorem for integration, 

j A p fd p x = A p f(x*) J d p x = A p f(x*) ■ vol(D), 



V 



V 



for some x* in D. But in the limit as the region of integration becomes 
infinitesimally small, it does not matter which x* we use in T> since / is 
continuous, so 



J A p f d p x A p f ■ vol(P) as vol(£>) 



0. 



v 



Thus 



A„f = lim 

vol(£>)->-0 



vol(P) 



dV 



d_l 
dn 



da 



Note that in this limit, the volume of T> approaches wr p_1 br. We can break up 
the integral in the numerator into one integral over the right bounding surface, 
one over the left bounding surface, and one over the lateral surface which we 
will denote by dT>' . Since in this limit the bounding surfaces are infinitesimally 
small, we can pull the integrands outside the integrals and rewrite the above 
equation as 



df 



A„ f = hm 



W\r+5M r+ Sr ) 



p-1 



9£\ 

dr I 



UJr 



p-1 



urP 1 8r 

which becomes, using the binomial expansion, 10 



+ lim 

vol(D)^0 



JdV' I dn 



da 



vol(X>) 



lim 

<5r->0 

or, 



Ul+Sr + (P ~ V 1 ^* 6r + °(^ 2 )] ~ wi-^ , contribution from 



A„ f = lim 



9f\ 


- ^1 


dr \r- 


- Sr dr \r 



Sr 



+ 



+ 



lateral surface, 



P~ 1 df_ 

r dr 



+ 



contribution from 
lateral surface. 



5 We introduce here our notation for a normal derivative: §i = V„ f ■ h. 
Here, 0(Sr ) denotes terms which are such that t °™ s remains finite as Sr 
particular, -> as Sr -> 0. 



10, 



(2.13) 



0. In 
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Notice that the contribution-from-the-lateral-surface term 



lim 

vol(Z>)->0 



1 



df 

vol(P) J dn 



da 



contains directional derivatives only in directions orthogonal to the radial di- 
rection. This term thus contains no derivatives with respect to r and only 
derivatives with respect to the angles (j), 9%, 02, ■ ■ ■ , p -2] we will indicate it by 

A S P~lf, 



where we determined the prefactor 1/r 2 through dimensional analysis. We 
call the operator Ag P -i the spherical Laplace operator in p — 1 dimensions. 
Equation (2.13) thus implies the following proposition. 



Proposition 2.5 



dr' 1 



+ 



V 



Id 1 



Or 



Now that we are comfortable working in W let us move on to briefly study 
the subject of orthogonal polynomials which will be very useful when we come 
to our main topic of discussion. 



Chapter 3 

Orthogonal Polynomials 



This chapter provides the reader with several useful facts from the theory of 
orthogonal polynomials that we will require in order to prove some of the 
theorems in our main discussion. In what follows, we will assume the reader 
has taken an introductory course in linear algebra and, in particular, is familiar 
with inner product spaces and the Gram-Schmidt orthogonalization procedure. 
If this is not the case, we refer him or her to [4]. 

3.1 Orthogonality and Expansions 

We will deal with a set V where the members are real polynomials of finite 
degrees in one variable defined on some interval 1 /. With vector addition 
defined as the usual addition of polynomials and scalar multiplication defined 
as ordinary multiplication by real numbers, it is easy to check that we have, 
so far, a well defined vector space. 

Let us put more structure on this space by adding an inner product. For 
any p,q £ V, define a function (•, -) w : V x V — > M. by 



where w, called a weight function, is some positive function defined on / such 



exists and is finite for all polynomials r G V . Notice that, since the product 
of any two polynomials is itself a polynomial, the above requirement ensures 

1 For now, we allow I to be finite or infinite and open, closed, or neither. 




(3.1) 



that 
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3. Orthogonal Polynomials 



the existence of {p, q) w for all polynomials p, q. Moreover, the function (3.1) is 
symmetric and linear in both arguments and, since w is positive, (r, r) w > 
for all polynomials r. Finally, since r 2 is continuous for all polynomials r, 

1 

For suppose r(xo) ^ for some xo € I. By continuity, r must not vanish 
in some neighborhood about xo, and thus the above integral will acquire a 
positive value. Therefore, V together with the function (3.1) is a well defined 
inner product space. We call (•, -) w the inner product with respect to the weight 
w. 

The reader should recall the following two inequalities that hold in any 
inner product space. We state them here for reference without any proofs 
which may be found in [4]. 

Proposition 3.1 (Cauchy-Schwarz Inequality) For any vectors p,q in an 
inner product space V, we have \ {p,q)\ < ||j?||[M|- 

Proposition 3.2 (Triangle Inequality) 2 For any vectors p, q in an inner 
product space V , we have \\p + q\\ < ||p|| + 

Now that we have an inner product space, we can speak of orthogonal poly- 
nomials. We call two polynomials p, q orthogonal with respect to the weight? 
w provided (p, q) w = 0. 

Suppose now that we have a basis B = {4> n }^L for V consisting of or- 
thogonal polynomials where <f> n is of degree n. We can easily see that such a 
basis exists by beginning with a set of monomials {x n }^L and, using the 
Gram-Schmidt process, to come up with the {4> n }^ =0 : We leave the first 
member x° = 1 of the basis as is. That is, <fio = 1. Then, from the sec- 
ond member x 1 = x, we subtract its component along the first member to get 
<fti = x — l^jy • 1, so that the second member of our new basis is orthogonal 
to the first. We continue in this fashion, at each step subtracting from x n its 
components along the first n — 1 orthogonalized basis members. This creates 
an orthogonal basis for V . 

Now since B is a basis, given any polynomial q G V of degree k we can 
write 4 , for some scalars c n , 

oo 

q = } ] Cn(j) n , (3.2) 

n=0 

2 This is also called the Minkowski inequality. 

3 Although the weight function is important, often we tend not to stress it and call the 
two polynomials simply orthogonal. We also use the simplified notation (■,■}. 

4 We do not need to worry about the convergence of the following series — we will see in 
a moment that it is actually a finite sum. 
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i.e., we can expand q in terms of the <p n . Since each polynomial <f> n in B has 
degree n, we can clearly write any q in terms of the set Bk = {4> n } n =o- Indeed, 
the space of all polynomials of degree k or less has dimension k + 1, and the 
linearly independent set Bk has k + 1 elements. Using the uniqueness of the 
expansion (3.2), we can see that 

c n = for all n > k = degq. (3-3) 

Since the coefficients are given by 

_ (q,(t>n) 

we have the following result. 

Proposition 3.3 Given any polynomial q and any set of orthogonal polyno- 
mials {0 n }^Lo> where 4> n has degree n, 

{4>n, q)w = J 4> n q wdx = 0, for all n > deg q. 
l 



Since we will come across expressions of the form (p,p) frequently, let us 
define the norm of p with respect to the weight w(x) as 

IIpIU = f V (p,p)w 

Keeping in mind the way we constructed our set of orthogonal polynomials 
{4> n }%Lo (requiring 0o to have degree — i.e., requiring <fo to be the constant 
polynomial 4>q = ao — and choosing each successive (j) n to be of degree n and 
orthogonal to the previous n members) the following result should not surprise 
the reader. 

Proposition 3.4 Any two orthogonal bases of polynomials {0 n }^L o and {i^ n } n ° =0 
for which deg 4> n = deg ij) n = n, must be the same up to some nonzero multi- 
plicative factors. 

Proof Since the {^nj^Lg constitute a basis, we can write any ip m as, using 
(3.3), 

m 
n=0 

for some constants c n given by 

Cn = II 1 Wl ^rni^n)- 
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By Proposition 3.3, c n = for all n ^ m. We are thus left with 

and we are done. ■ 

3.2 The Recurrence Formula 

We can find a useful recursive relationship between any three consecutive 
orthogonal polynomials n +i, <j> n , (j> n -i- I n what follows, we will let £ n 
denote the coefficients of the x n term and the x n_1 term in <p n , respectively. 

Proposition 3.5 Given any set of orthogonal polynomials B = {0 n }^_ O ; 
where (j) n has degree n, 

4> n+ i - (A n x + B n )<t> n + C n <f) n -i = 0, (3.4) 



k n+l a ( Zn+l £n\ n A n \\(j) n \\ 2 



where 



k ' n \h ,i h ' A , Wch Jl 2 ' 

»>n V^n+l ™n/ \\Yn- 1\\ 

and Co = 0. 

Proof From the formula given for A n in (3.5), we see that (f> n +i — A n x4> n must 
be a polynomial of degree n since the leading term of ^> n +i is canceled. Since 
the set B is a basis, we can thus write 

n 

(j) n+ l - A n x<j) n = y^7j(/>j, 

3=0 

using (3.3). Taking the inner product of each side with c/y, where < f < n, 
we find 

7i ' = ||^.,||2 (i^n+lAf) ~ A n (x<p n , (j)j<)) , 

using the orthogonality of B and the linearity of the inner product. From the 
definition of the inner product (3.1), we see that the above equation is the 
same as 

7j ' = pV|j2 - A n{<t>n, X(f>j')) ■ 

Since f ^ n + 1, the first inner product in the above equation vanishes. Since 
deg(xi^y) = deg(</y) + 1 = j' + 1, we can use Proposition 3.3 to determine 
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that the second inner product above is zero unless f = n — 1 or j' = n. If we 
rename -y n = B n and 7 n _i = — C n , we have 

4>n+l ~ {A n X + B n )(f) n + C n (f>n~l = 0, 

as in (3.4), with 

—A A 

B n = 1, 7 i?o (4>n,x4> n ), C n = " (0n> E0n-l), (3-6) 

and it remains to compute the coefficients B n ,C n . We will start with C n , 
rewriting x(j) n -i = x(k n -\x n ~ 1 + • • • ) as 

k -i 

k n -\x n + lower order terms = — — (k n x n + lower order terms) 

fen 

' + lower order terms) . 



k n 



Since the lower order terms 'die' in the inner product with <p n , (3.6) implies 



C n 



A n fen— 1 ||^ || 2 A n || (f> n || 2 

f'n-ill 2 K n A n _i ||0 n _i|| 2t 



as required. We calculate B n in a similar way, rewriting x<p n = k n x n+1 + 

^■f^X I * 3iS 



X0 n 



a; n + L.O.T. 



kn+l&n 
k 

k n x n + L.O.T. 
6 n + L.O.T. 



where L.O.T. is an acronym which stands for "lower order terms". Only the 
(f) n term will survive the inner product with 4> n , and (3.6) implies 

n ~A n f &n 4+1 1 ii i n2 a ( ^n 
B n = -n— 7T7T ~, \\0 n \\ = A, 



kn 


k n+ lX n+1 + 


i n+ ix n 


k n +i 


kn 


+ {Vn 




k n +l 


k n +i 


kn 






k n +l 


k n +i 



^n|| 2 \k n fen+1 7 \fen+l k, 

and we are done. ■ 

We see that if we know the coefficients A n , B n , C n then once we have chosen 
two consecutive members of our set of orthogonal polynomials, the rest are 
determined. 
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3.3 The Rodrigues Formula 

Until now, we spoke of the arbitrary interval /. For the remainder of our 
discussion, we will be most concerned with the interval [—1,1]. 

Consider the functions ifi n (%) defined on [—1, 1] for n = 0, 1, . . . as 

*» w = ^(^)"N*)(i-* 2 )"]- (3 ' 7) 

For arbitrary choice of the weight function w, we cannot say whether these 
functions are polynomials nor whether they are orthogonal with respect to w. 
In an effort to force the ip n to be polynomials, we can start by looking at tpi 
and requiring it to be a polynomial of degree 1, 

, , . w'(x) 9 . „ roq 

ifj^x) = - x 2 ) - 2x = ax + b. 

w(x) 

This creates a differential equation, which we can write as 

w' _ (a + 2)x + b 
w (1 + x)(l — x) 

Integrating, 

In w = f 7^ — ) — — — r dx + const. 
J il + x){l-x) 

Decomposing the integrand into partial fractions to compute the integral, we 
find 

a + b + 2 i . , b-a-2, 
mw = in (1 — x) H In (1 + x) + const., 

or, 

w(x) = const, x (1 — x) a (l + x) 13 , 

where a, (3 are the coefficients from the previous expression. Clearly, the ar- 
bitrary constant in the above expression cancels when w{x) is inserted into 
(3.7), so we will take it to be unity. We will also require a, ft > — 1 so that we 
can integrate any polynomial with respect to this weight. For our purposes, 
we will see that these restrictions on w{x) are all we need. 

Proposition 3.6 Let w(x) = (1 - x) a (l + x) 13 with a, (3 > -1. Then (3.7) 
defines a set of polynomials {ipn}%Lo wnere degip n = n. 

In what follows, we will use (k)i to denote the falling factorial, 
{k)t = f k(k - 1) • • • (k - I + 1) for each t G N, 
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and 



n\ def nl 



k) k\(n-k)\ 



Proof We can see this using the Leibnitz rule for differentiating products, 



*X(f )=T( n ) dkf dn ~ kg 

dx J [ ' ^\k)dx k dx n - k ' 
7 fc=o v 7 



which can be proved easily by induction. 

With the weight w(x) given previously, (3.7) becomes 



t \ n 

M*) = (1 - »)" a (l + ^ ( "H (1 " ^) a+n (! + z) 



(i.r 



"" ' — )(l-^ n+a 



(1 + z)^ ( ~ ) (1 + z) n+ 



EC) d--) 

fc=0 

(-l) k {n + a) k {l-x) n - k (n + /3) n _ fc (l + x) 



n— fc 



fc=0 



We can see that ipn is a polynomial. Let us examine the leading term in tpn- 
We see from the last line in the above equation that 



k=0 ^ 7 



(n + a) k (-x) n - k (n + P) n - k (xY + L.O.T. 



= (-l)V^(^(n + a) fc (n + /3) n _ A; + L.O.T. 
fc=o ^ 7 

Since the sum in the last line of the above equation is strictly positive, we see 
that the x n term in tp n has a nonvanishing coefficient; i.e., the degree of ip n is 
precisely n. ■ 

Proposition 3.7 Let w(x) = (1 — x) a {\ + Xf , for a, (3 > —1. Then the ip n 
defined in (3.7) satisfy 



ip n (x)x k w(x) dx = 0, for < k < n. 
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Proof Let J be the above integral. Inserting (3.7) for ip n and integrating by 
parts n times, 



J 



d 

dx 



[w{x)(l - x 2 ) n ] dx 



dx 



n-l 



[w(x)(l - x 2 ) n ] 



kx 



k-l 



dx 



n-l 



[w(x)(l -x 2 ) n ] dx 



Y^(-l) e (k) e x 



k-l 



1=0 



dx 



n-l-l 



[w(x)(l - x 2 ) n ] 



Upon inserting the assumed form of w(x), this becomes 



J = Y / (-i) e (k) 



f;X 



k-l 



1=0 



d_ 
dx 



n-i-l 



{l-x) n+a {l + x) 



n+/3 



In each term of the above sum, the factors (1 — x) n+a and (1 + x) n+l3 keep 
positive exponents after being operated upon by (^)" 1 so that J vanishes 
after we evaluate the sum at —1 and 1. ■ 

Since each ip n has degree n and is orthogonal to all polynomials with degree 
less than n, we have the following result. 

Corollary 3.8 The formula (3.7) defines an orthogonal set of polynomials 
{V'nj^Lo known as the Jacobi polynomials with degip n = n for each n. 

We have thus found that (3.7), called the Rodrigues formula, with an ac- 
ceptable weight function w(x), defines a set of orthogonal polynomials. We 
are free to multiply each polynomial ip n by an arbitrary constant without up- 
setting this fact. The different classical orthogonal polynomials can be defined 
using the Rodrigues formula by adjusting this constant and the exponents a, j3 
in w{x). 

Now, suppose we discover that a set of polynomials {</> n }^l such that 
deg 4> n — n is orthogonal with respect to the weight w(x). Can we conclude 
that these polynomials satisfy the Rodrigues formula (3.7)? Comparing Corol- 
lary 3.8 and Proposition 3.4, we can answer "yes." That is, we can define the 
(f) n by the Rodrigues formula with suitable multiplicative constants, 



<Pn(x) 



d 

w(x) \dx 



[w(x)(l 



3.4. Approximations by Polynomials 



37 



3.4 Approximations by Polynomials 

The main result in this section will be the Weierstrass approximation theorem, 
which will allow us to use polynomials to approximate continuous functions 
on closed intervals. First, we will introduce some ideas that we will use to 
prove this important theorem. 

Given a function / : [0, 1] — > M, we will use B n (x; f) to denote a Bernstein 
polynomial, defined on the closed interval [0, 1] as 



^ /) = £(I)/ 

fc=0 v J 



- ) x k {l-x) n -< 



n 



for n = 1,2,.... Since / (^) is just a constant, these are clearly polynomials. 
We will compute three special cases of the Bernstein polynomials. 

Lemma 3.9 For every natural number n, 

B n (x;l) = l, (3.8) 
B n (x; x) = x, (3.9) 

B n {x-x 2 ) = x 2 + X{1 ~ X \ (3.10) 
n 

Proof By the binomial theorem, 

k=0 ^ ' 

Substituting y = 1 — x gives 

1 = E (D - x ^~ k = B n(x; 1), (3.12) 

k=0 W 

proving (3.8). Differentiating (3.11) with respect to x, 



fc=0 ^ ' 

and multiplying by x/n, 

x(x + y)"" 1 = £ (jf) (I) xk V n ~ k - ( 3 - 13 ) 
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Substituting y = 1 — x gives 



x = 

k=o 



proving (3.9). Differentiating (3.13) with respect to x, 

(x + y)"" 1 + (n - !)*(* + yf" 2 = £ Q (£) 



2^ 



fc=0 

and multiplying by x/n, 



2 

k„.n—k 



(^)(-^£(:)00 * 



fc=0 

Substituting y = 1 — x gives 

* 2 + = t (I) (£) 2 " = x 2 ), 

proving (3.10). ■ 

We are now well equipped to prove the following result. 

Theorem 3.10 (Weierstrass Approximation Theorem) Let the function 
f : [a,b] — > M. be continuous, and let e > 0. Then there exists a polynomial 
p(x) such that \f{x) — p(x)\ < e for all x in the interval [a,b]. 

We will show that the Bernstein polynomials, for sufficiently large n, work 
to approximate the function / to the required accuracy. 

Proof Let e > 0. First, we will redefine the function / so that we can work 
in the interval [0,1]. Using the linear transformation 

T(x) = (b — a)x + a, 

we set 

/ = / o T : [0, 1] R. 

If we can find a polynomial p to adequately approximate /, we can use the 
polynomial poT^ 1 : [a, b] — > R as our required approximation of /. From here 
on, we will write / instead of / for simplicity. 
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Since / is continuous on the closed and bounded interval [0, 1], we know 
that / is bounded and uniformly continuous. By the definition of boundedness, 
there exists an M such that 

\f(x)\ < M, for every x G [0, 1], 



which implies 
!/(*)-/(*) 



77 



< 2M, for all x £ [0, 1] and < k < n. 

(3.14) 



By the definition of uniform continuity, there exists a 5 > such that 



f(x) - f 



n 



< -, whenever \x 1 < 5. 

2 n 



(3.15) 



Now, let 

E=\f(x)-B n (x;f)\, 

and we will estimate this error from above. Using (3.8) and the triangle 
inequality, 



H'w-£(r>(O x ' (i 

k=0 v 7 v 7 

n 



\n—k 



k=0 



f(x) ~ f 



n—k 



x\l-x) 



x k (l-x) n - k . 



Splitting up the sum, 

E 

\x-±\<5 

1 n 1 

+ E 

i*-£i>* 

Using (3.14) and (3.15), 



n—k 



f{ X )-f[-) x\l-x) 



n—k 



E <\ E (^(i-xr*+2M e ftW 

I™ I ^ r I ™ I \ £ 



\n—k 



\x-*-\<8 

1 r). 1 



1 T). 1 



" 2E 1 • 

fc=0 



A; 



d-x)- + ^E(;) 

fc=0 v 7 



< J x fc (l-x) n - fc , 

n , 
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or 

E < \ + ~W [^"(^ 1 ) ~ 2xB n(x; x) + B n (x; x 2 )] 
Using Lemma 3.9, 

e 2M / o n 9 x(l - x) 

^ < o + -«r \ x 2x + x + 

2 \ n 

Since the function x(l — x) has a maximum value of 1/4, we can write 

„ e 2M p „ M 

2 4n£ 2 e<5 2 

Thus, we see that the Bernstein polynomials do the job for large enough n, 
completing the proof. ■ 

Now that we know polynomials are a good tool for approximating contin- 
uous functions, we ask how to find the best such approximation. This is a 
result from linear algebra. For a function /, we define the best approximation 
p to be the one that minimizes the norm of the error, i.e., for which ||/ — p\\ 
is smallest. 

Proposition 3.11 Let f be a function, and let {^fcl^Lo be an orthogonal set 
of polynomials with deg(p n = n. The polynomial 

Pn = y^a k <p k , where a k = ' , (3.16) 
fro H^ll 

is the unique polynomial of degree n that best approximates f , i.e., that mini- 
mizes ||/ — p n \ . 

Proof Let's choose a different polynomial p n of degree n. We will write p n as 

n 

Pn=Pn + Qn, where q n = b k (p k , 

k=0 

for some constants b k not all zero. Then, we can write ||/ — p n \\ 2 as 

11/ ~Pn ~ Qn\\ 2 = ||/ -Pnf + hnf ~ 2(f -p n ,q n ). 

But 

n 

(f ~ Pn, Qn) = ^2 b k ((/, 4>k) - (Pn, 4>k)) = 0, 
k=0 

since each term in parentheses vanishes, as we can see from (3.16). Thus, for 
any p n , 

ll/-Pn|| 2 = ||/-Pn|| 2 + |kn|| 2 > ||/-Pn|| 2 , 

the error is greater than that of p n . ■ 
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This fact is really just a special case of a more general result from the study 
of inner product spaces in linear algebra. We will state the more general fact 
next. The proof is exactly the same. 

Proposition 3.12 Let f be a vector in an inner product space, and let {^fcjfcLo 
be an orthogonal set of basis vectors. The vector 



Pn = '^2a k (t>k, where a k 



is the unique linear combination of the first n + 1 basis vectors that best ap- 
proximates f, i.e., that minimizes \\f — p n \\. 



3.5 Hilbert Space and Completeness 

Later in the main discussion, we will be interested in expanding an arbitrary 
function / in an infinite series of the spherical harmonic functions. In this 
section, we will address the question of when such an expansion is possible. 
Since we are concerned primarily with spherical harmonics in an arbitrary 
number of dimensions, we will keep the discussion general, making no reference 
here to the number of variables on which the functions depend. In what follows 
we will thus let x denote a vector in W, V denote a subset of W, and dQ denote 
the differential volume element of T>. 

We will restrict our attention to the space of square-integrable functions 
with domain T> with respect to the weight w(x), i.e., those functions / : T> — > M 
for which the integral 

J f(x) 2 w(x)dn 

exists and is finite. When endowed with the usual operations of function 
addition and multiplication by real numbers and paired with the inner product 

(f,g) = J f(x)g(x)w(x)dn, 

the set of such functions becomes an inner product space. Such a space has 
an induced norm 

11/11 = y/UJ). 
We will now work up to a definition of a Hilbert space. 

Definition In a normed space, a sequence {x n }^ =0 is called Cauchy provided 
that for any e > there exists an N such that 

\\x n — x m \\ < e for all n,m> N. 
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A Cauchy sequence is always convergent. However, it is possible that its 
limit does not belong to the same space. 

Definition A normed space is complete provided every Cauchy sequence con- 
verges to an element in the space. 

Definition A complete inner product space is called a Hilbert space. 

We will state the following fact without proof. The interested reader may 
consult [9] for the proof. 

Theorem 3.13 The inner product space of square-integrable functions defined 
above is a Hilbert space. 

We can restate the question that opened this section as follows. Given a 
Hilbert space % and an orthonormal set {^„}^L C % (which can be obtained 
from any orthogonal set by dividing each member by its norm so ||</>n|| = 1), 
when can we write any arbitrary member / S H as a linear combination 5 of 
the (/>„? 

Definition An orthonormal set {(f> n }^Lo f= 7~L is called complete provided that 
for each / 6 T~L, there exist scalars ci, C2, . . ., such that 

n 

lim ||/- J^cM =0. (3.17) 

k=0 

We know from Proposition 3.12 that out of all linear combinations, the 
combination 

n 

Pn = ^{f, 4>k)(l>k 
k=0 

minimizes ||/ — p n \\. Notice that we have now normalized the vectors 4>k to 
unity: \\4>k\\ = 1- Thus, if there exist scalars ct for which the norm in (3.17) 
converges to zero, then certainly 

n 

lim ||/- V(/,^)0 fc || =0, (3.18) 

k=0 

since 

n n 

11/ - ^2{f,<j>k)<f>k\\ < 11/ - 5^Cfc0fe||, for every n. 

k=0 k=0 

A linear combination is usually assumed to contain a finite number of terms. We will 
not need this restriction for our purposes. Here, linear combinations can be finite or infinite 
sums. 
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9fc/<?>fc 



Let us rewrite (3.18) by computing 

n 

/-£</« 

which is the same as 

n 

2(f, 4> k )<p k ) + ( fc }<^, ]T</, 



fc=0 



fc=0 fc=o 
Using the linearity of the inner product, this becomes 



k,9t 



k=0 k,i=0 
By the orthonormality of the n , we have 



<>k)(Pk 



k=0 



f ~ £</», 
k=0 

We note that since the norm-squared is non- negative, then 

n 

-^(/,^) 2 >o, 



fc=0 



or, 



£</,< 

k=0 



lk?< 



for every n. 



(3.19) 



(3.20) 



This is known as BesseVs inequality and tells us that the sum Ylk^oif^k} 2 
converges. If the <f>k form a complete set, ||/— Y^=oifn, ^fc)^fe|| 2 must converge 
to zero as n — > oo. In this case, from equation (3.19), it follows that Bessel's 
inequality becomes an equality: 

oo 
fc=0 

known as Parseval's equality. We have thus arrived at the following conclusion. 

Proposition 3.14 An orthonormal set {</>n}^Lo is complete if and only 
if Parseval's equality holds for each f G %. 



We give one more definition. 
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Definition A set {<j> n }'%Lo is closed provided 

(/> 0n) = for every n implies / = 0. 

Now we are ready to prove the main result of this section. 

Theorem 3.15 The orthonormal set {(/> n }^L Q7i is complete if and only if 
it is closed. 

Proof (Sufficient) Let {4> n }^ =0 C H constitute a closed orthonormal set and 
/ G H. We will show that 



lim ||/f-^(/,^) 2 =0, 



fc=0 



since this guarantees completeness by Proposition 3.14. 

Define g n = f - Yl=o{fAk)4>k, and notice that {g n }%L is a Cauchy 
sequence since, if n > m, 



\\9n ~ 9m\\ 2 = || ^2 (f,(f>k)(t>k 
k=m+l 
n 

= E U^k){fAi>){4>kAi) 

k,£=m+l 
n 

fc=m+l 

which can be made arbitrarily small for large enough n, m since the series 
Yl'kLoif'^k) 2 converges, by (3.20). Since H is complete by definition, this 
Cauchy sequence must converge to some g S H, i.e., 

lim \\g n -g\\ =0. (3.21) 

n— >oo 

Now fix j and take n £ No such that n > j. By the Cauchy-Schwartz inequality 

\(9,4>j}\ = \(9n - 9,4>j}\ < \\9n ~ g\\\\(t>j\\ = \\9n ~ 9\\- 

Since this holds for any n > j, we can use (3.21) to conclude 

\{g,<t>j)\ < lim hn - g\\ = o, 

which implies that for arbitrary j, (g,(f>j) = 0. Since we assumed the cftj 
constitute a closed set, we have g = 0. Therefore, 

lim ( Il/H 2 - V(/,0 fc ) 2 ] = lim \\g n \\ 2 = lim \\g n - g\\ 2 = 0. 
\ k=0 / 
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(Necessary) Suppose {4> n }'^ 3 =0 is complete but not closed. Then there exists 
a function / 7^ such that (/, cp n ) = for every n. Then 

Hm (||/|| 2 -£</,0 fc > 2 ) =||/|| VO. 
V fc=o / 

So, by Proposition 3.14, {(pn}^^ is not complete — a contradiction. This 
proves that if an orthonormal set in % is closed, then it is complete. ■ 

In Chapter 4, we will use this result to show that the set of spherical 
harmonics form a complete set in the space of square-integrable functions by 
showing that it is closed. 



Chapter 4 

Spherical Harmonics in p 
Dimensions 



We will begin by developing some facts about a special kind of polynomials. 
We will then define a spherical harmonic to be one of these polynomials with a 
restricted domain, as hinted at in Section 1.1. After discussing some properties 
of spherical harmonics, we will introduce the Legendre polynomials. And once 
we have produced a considerable number of results, we will move on to an 
application of the material developed to boundary value problems. 



4.1 Harmonic Homogeneous Polynomials 

Definition A polynomial H n (xi,x 2 , . . . ,x p ) is homogeneous of degree n in 
the p variables x\, X2, ■ ■ ■ , x p provided 

H n {tx 1 ,tx 2 , . . -,tx p ) = t n H n (x 1 ,x 2 , ...,x p ) . 

In the definition of a homogeneous polynomial, let's set U{ = txi, for all i 
and differentiate the defining equation with respect to t: 

EdH n (u 1 ,u 2 , . . . ,u p ) dui ! . 
o ^r= nt H n {x 1 ,x 2 ,...,x p ), 
oui dt 

i=i 

or 

J^dH n ( Ul ,u 2 ,...,u p ) x 

} J q- ^ Xi = nt H n (x 1 ,x 2 ,...,x p ). 

i=l U% 

Finally, we set t = 1 to find the following functional equation satisfied by the 
homogeneous polynomial, 

EdH n (xi,x 2 , ...,x p ) 
Xi = nH n (xi,x 2 ,. . . ,x p ) , (4.1) 
dxi 

i=i 
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known as Euler's equation. 

The following calculation will be useful in the counting of linearly inde- 
pendent homogeneous polynomials. 



Lemma 4.1 For < |r| < 1, 



( l_ r)P j!(p_i)! 



Proof We will use a counting trick to prove this. Since < \r\ < 1, we can 
write (1 — r)~ p in terms of a product of geometric series, i.e., 

(-)--(^)(r^)-(i^) 

V v ' 

p times 

(oo \ / oo \ /oo \ 
E r " E r " ••• E r " • 
n=0 / \n=0 / \n=0 / 
v v ' 

p times 

We will compute the Cauchy product of the p infinite series. The result will 
be an infinite series including a constant term and all positive integer powers 
of r, 

oo 

(l-r)-" = J>ri. 

3=0 

It remains to compute the coefficients Cj. To determine each Cj, we must 
compute how many r J 's are produced in the Cauchy product, i.e., in how 
many different ways we produce an r° in the multiplication. 

Note that this computation is equivalent to asking, "In how many differ- 
ent ways can we place j indistinguishable balls into p boxes?" The p boxes 
correspond to the p series in the product, and choosing to place k < j balls 
into a certain box corresponds to choosing the r k term in that series when 
computing a term of the Cauchy product. 

Let us use a diagram to assist in this calculation. We will use a vertical 
line to denote a division between two boxes and a dot to denote a ball. Each 
configuration of the j balls in p boxes can thus be represented by a string of j 
dots and p — 1 lines (since p boxes require only p — 1 divisions). For example, 



represents one configuration of j = 8 balls in p = 8 boxes. 
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Now, the number of ways to arrange j indistinguishable balls in p boxes 
is the same as the number of distinct arrangements of j dots and p—1 lines. 
This is given by u p — 1 + j choose j" i.e., 

p-l+j\ _ (p- 1 + j)l 



j J j\(p-l)\ 



and the lemma is proved. 



Proposition 4.2 // K (p, n) denotes the number of linearly independent ho- 
mogeneous polynomials of degree n in p variables, then 

K ^ n > n\(p-l)\ ' 



We will give two proofs of this claim. The first uses a recursive relation 
obeyed by K(p,n), while the second employs another counting trick. 

Proof 1 Let H n (xi,X2, • • • , x p ) be a homogeneous polynomial of degree n in 
its p variables. Notice that H n is a polynomial in x p of degree at most n. For 
if H n contained a power of x p greater than n, the polynomial could not be 
homogenous of degree n, since H n {. . . , tx p ) would contain a power of t greater 
than n. Thus we can write, 

n 

, X2 , • • • , 3?p) — ^ X p fl n —j (xx , X2, ■ ■ ■ , Xp—l 

), (4-3) 

3=0 

where the h n -j are polynomials. Moreover, notice that the h n -j must be 
homogeneous of degree n — j in their p — 1 variables. Indeed, using the homo- 
geneity of H n , 

n n 

t ^ x p h n —j (xi , . . . , — H n (tx\ , . . . , iXp) — ^ j (t%p)^ hn—j itx\ tx p —\) , 

3=0 3=0 

which implies that 

n 

[t n ~ j h n -j(xi, . . . ,x p -i) - hn-j(txi, . . . ,txp-{)\ = 0. 

3=0 

The expression in brackets must vanish by the linear independence of the 
xp. Each h n -j can be written in terms of a basis of K(p— 1, n—j) homogeneous 
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polynomials of degree n — j in p — 1 variables, and thus H n can be written in 
terms of a basis of 

n n 

K(p,n) = ^K(p- l,n- j) = ^K(p- 

3=0 3=0 

linearly independent elements. We have found a recursive relation that the 
K(p,n) must satisfy. Now for some < |r| < 1, let 

oo 

G(p) = ^Tr n K(p,n). (4.4) 



Then, 



or 



so 



n=0 



oo n oo oo 

g( P ) = j2 rn T, K (p- 1 >ri=H K (p -ME r "> 

n=0 J=0 J=0 



oo oo 1 oo 

g(p) = y, K (p - 1. jV E r " = r;E ^(p - x ' ■?')' 

j=0 n=0 j=0 



1 - r 

Using an inductive argument, we can show that 

G(l) 



G(p) 



(l_ r )P-i • 



rt in one variable can be written as ex™, we see that 



By noticing that K(l,n) = 1, since every homogeneous polynomial of degree 

ritten as ex™, we see 

oo 

G(l)= j>» = — 



3=0 

Thus, 

G{p) = (l-r)-v = ± {p + n - l K \ 
xr> y ' ^ n\(p-l)\ 

n=0 ^ ' 

by the above lemma. Comparing this to (4.4) then proves the theorem. ■ 

Proof 2 Using reasoning similar to that used in the above proof, we see that 
every monomial (i.e., product of variables) in an n-th degree homogeneous 
polynomial must be of degree n. For example, a fourth degree homogeneous 
polynomial in the variables x, y can have an x 2 y 2 but not an x 2 y term. 
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To uniquely determine such a polynomial, we must give the coefficient of 
every possible n-th degree monomial. So to find K(p,n), we must find out 
how many possible n-th degree monomials there are in p variables. But this 
is equivalent to asking, "In how many ways can we place n indistinguishable 
balls into p boxes?" The p boxes represent the p variables from which we can 
choose, and placing j < n balls into a certain box corresponds to choosing to 
raise that variable to the j power. For example, using the same notation as 
in the proof of Lemma 4.1, the string 



could represent the w 2 x 1 y 1 z° term in a fourth-degree homogeneous polyno- 
mial. As in the above lemma, the number of such arrangements is 

(n+p-l\ (n + p-iy 
K(p, n) - 1 1 



n J n\{p — 1)! 
proving the theorem. ■ 

Definition A polynomial q(xi,X2, ■ ■ ■ ,x p ) is harmonic provided 

A p q = 0. (4.5) 

Equation (4.5) is called the Laplace equation. 

The following property of combinations will be used to prove the theorem 
that follows it. 

Lemma 4.3 Ifk,£<E N, then 

fk\ _ k fk-l 
\£J ~ ~£\£- 1 

Proof Just compute: 

kfk-l\ k (k-l)l k\ (k 



\£-lJ £ (£-l)\(k-£)\ £\{k-()\ \£ 



as required. 



Theorem 4.4 // N(p, n) denotes the number of linearly independent homo- 
geneous harmonic polynomials of degree n in p variables, then 

at/ \ 2n+p-2/n + p-3 
Jy (p,n) = - 



n \ n — 1 
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Proof Let H n be a homogeneous harmonic polynomial of degree n in p vari- 
ables. As in (4.3), we write 

n 

H n (xi,X 2 , ...,X p ) = y~]x 3 p h n -j(xi,X2, ■ ■ .,Xp-l), 
3=0 

and we operate on H n with the Laplace operator. Since H n is harmonic, 

/ d 2 \ n n 

= A p H n = + A p _iJ H n = y~]j(j - l)x J p ~ 2 h n -j +^x p A p - 1 h n -j, 



so 



= J24 tO' + 2 )0' + l ) h n-j-2 + Ap-lhn-j] . 

where we define h—\ = /i_2 = 0. Since the xp are linearly independent, each 
of the coefficients must vanish, 

2/ln-2 + \-lK = 0, 
6K-3 + Ap-lhn-l = 0, 

: (4-6) 
n(n - l)h + A p _i/i 2 = 0, 
Ap_x/ii = 0, 
Ap_i/io = 0. 

We have found a recursive relationship that the h n ~ must obey. Thus, 
choosing /i n and /i n -i determines the rest of the h n -j. By Theorem 4.2, 
h n can be written in terms of K{p — l,n) basis polynomials and /i n _i can 
be written in terms of K (p — 1, n — 1) basis polynomials. Thus we must give 
K(p— 1, n) + — 1j n ~ 1) coefficients to determine /i n , /i n _i which thereby 
determine H n . Therefore, 

N(p,n)=K(p-l,n) + K(p-l,n-l)= (^ + ™~ 2 ) + ( p + n ~ 3 

By the above lemma, we can write this as 



n J \ n — 1 



p + n — 2/p + n — 3\ fp + n — 3\ 2n + p — 2/ra + p — 3 N 

iv(p, n) = £ r 1 + " ] ' 



n V n — 1 / \ n — 1 / n V n — 1 



and thus the theorem is proved. 



4.2. Spherical Harmonics and Orthogonality 



53 



In what follows, x will always denote the vector (x±, X2, ■ ■ • , x p ), r or 

|x| will denote its norm \J x\ + x\ + • • • + x^ , and £ will denote the vector 

^2> • • • j Cp) having unit norm. Keep in mind that x,£ represent vectors 
while the denote their components. Also, H n will always denote a har- 

monic homogeneous polynomial. Lastly, in the remainder of this chapter, we 
will often suppress explicit reference to the number of dimensions in which 
we work. It should be assumed that we deal with p dimensions unless stated 
otherwise. 



4.2 Spherical Harmonics and Orthogonality 

We are now ready to introduce the spherical harmonics. First, notice that 

H n (x) = H n (rO=r n H n (0. (4.7) 

In p-dimensional spherical coordinates, the radial dependence and the angular 
dependence of the functions H n can be separated. 

Definition A spherical harmonic of degree n, denoted Y n (£), is a harmonic 
homogeneous polynomial of degree n in p variables restricted to the unit (p— 1)- 
sphere. In other words, Y n is the map 

Y n : S p - 1 -> R, given by Y n (£) = H n {0 for every £ G S^ 1 

for some harmonic homogeneous polynomial H n . We can write Y n = H n \ SP ~i. 

In Chapter 1, we introduced functions Y (which we claimed were spherical 
harmonics) that turned out to be eigenfunctions of the angular part of the 
Laplace operator. Using Proposition 2.5, we can show that spherical harmonics 
are indeed eigenfunctions of Ag P -i. 

Proposition 4.5 

A SP -iY n = n{2 - p - n)Y n . (4.8) 

Proof Let H n be a harmonic homogeneous polynomial and Y n its associated 
spherical harmonic. As in (4.7), we have H n (x) = r n Y n {^). Then, using 
Proposition 2.5, 

= A p (r n Y n ) = n(n - l)r n ~ 2 Y n + V — nr n ~ x Y n + \ r n A SP - X Y n . 
Rearranging, 

r n - 2 [A SP -.Y n + n{n+p- 2)Y n ] = 0, 

which implies 

A SP -iY n = n(2-p - n)Y n , 
as sought. ■ 
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Remark In three dimensions (4.8) becomes 

„ d \ 1 d 2 



1 d 
sm6d0 \ 8m "d0 



+ 



sin 



Y, 



+ l)Y e . 



As promised in Subsection 1.2, this allows us to show that three-dimensional 
spherical harmonics carry a definite amount of quantum mechanical angular 
momentum. Referring to (1.13), 



L 2 Y e = h 2 £(£ + l)Y £ . 



Remark Since the spherical harmonic Y n (£) is defined as the restriction of 
some H n {x) to the unit sphere, Y n {£) must also be homogeneous. However, 
Y n (t£) is not always defined since the domain of the spherical harmonic is 
In fact, it is only defined for \t£\ = 1, and since |£| = 1, we see that 
we can only write Y n (t£) for t = ±1. The case t = 1 is trivial, but the case 
t = — 1 tells us the parity of Y n (£). Since Y n {— £) = (— l) n y n (£), we see that 
the transformation £ i — > — £ sends Y n i — — y„ if n is odd and leaves Y n 
invariant if n is even. 

We now come to the main result of this section, which we hinted at in 
Section 1.1. 

Theorem 4.6 Let Y n (£),Y m (£) be two spherical harmonics. Then 
J Y n (OY m (Odn p ^ = 0, ifn^m. 

5p-i 

That is, spherical harmonics of different degrees are orthogonal over the sphere. 



Proof Let us perform the computation using the harmonic homogeneous poly- 
nomials H n and H m where Y n = H n \g P -i and Y m = H m \ SP -i. 

Now, we start with the divergence theorem in p dimensions for a vector A, 



JV p -A{x)a?x= J AiO-ZdQj,-!, 

Bp 5p-i 
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apply it twice for the vectors H n VH m and H m VH n , and subtract the results 1 : 
J V p -[H n (x)V p H m (x) - H m (x)V p H n (x)] d p x 

BP 

= J [H n (OV p H m (0-H m (OV p H n (Z)]-Zdn p ^, 



Sp- 

or 



0= J [H n (^)V p H m (0-H m (OV p H n (0]-Cdn p ^, (4.9) 



S»-i 

where we used the property 



V p • (H n V p H m ) — V p H n • V p H m + H n A p H m , 



and the fact A p H m = A p H n = 0. 



Euler's equation (4.1) for a homogeneous polynomial, 



v 



— — Cj = 
i=l ^ 

can be written in the form 

V p ff„(0-£ = nfln(0- 
With the help of this result, equation (4.9) takes the form 

(m-n) y H n ^)H m (Odn p -x = 0. 

sp- 1 

But the integral is carried out over S' p_1 where y n = H n and y m = H m . This 
equation is thus equivalent to 

(m-n) j Y n {()Y m (Odn p - 1 = 0. 

sp- 1 

By hypothesis, n ^ m; therefore Y n , Y m are orthogonal over the sphere. ■ 



1 Incidentally, we point out that the identity thus obtained for any two functions /, g 
J (fA p g - gApf) (fx = J (fV p g - gV p f) ■ f dn p _i, 



sp- 1 



is known as Green's theorem. We may also observe that V p / ■ £ is the directional derivative 
of / along £ and write 

J(fA p g-gA p f)d?x= J ff|| - fl |0 dfi p _r. 
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Given a set of N(p, n) linearly independent spherical harmonics of degree 
n, we can use the Gram-Schmidt orthonormalization procedure to produce an 
orthonormal set of spherical harmonics, i.e., a set 

{YnAO}?Ji' n) with J Y n>i (Z)Y n j(£)<m p - 1 =6 ijl (4.10) 

where 

_ J 1 if % = j, 
" I iff /j, 

is the Kronecker delta. For the remainder of this chapter, unless indicated oth- 
erwise, we will let Y n ^{^) denote an n-th degree spherical harmonic belonging 
to an orthonormal set of N(p,n) such functions, as in (4.10). 

In what follows, let R be an orthogonal matrix that acts on £ as a rotation 
of coordinates. Notice that since the integration is taken over the entire sphere 
in (4.10), the orthonormal set of spherical harmonics remains orthonormal in 
a rotated coordinate frame. That is, 

J Y nA (Rt)Y nJ (RZ)cm p _ 1 = S ij . (4.11) 
Sp-i 

Proposition 4.7 IfY n (£) is a spherical harmonic of degree n, then Y^(£) = 
Y n (R£) is also a spherical harmonic of degree n, for any rotation matrix R. 

Proof Let Y n (£) be a spherical harmonic of degree n. Then there exists a 
harmonic homogeneous polynomial H n (x) of degree n such that Y n = H n \ SP -i . 
Denote H' n (x) = H n (Rx). We claim that Y^ = H' n \ SP -\. To see this, first 
notice that H' n (x) is a polynomial in aci, X2, ■ ■ ■ , x p . Indeed, H' n {x) = H n (Rx) is 
a linear combination of powers of the Y^j=i ^ij x j an< ^ thus a linear combination 
of powers of the Xj. Next, notice that H' n is homogeneous of degree n, 

H' n (tx) = H n {tRx) = t n H n {Rx) = t n H' n (x). 

Finally, notice that H' n is harmonic, by Proposition 2.1. Restricting H' n to the 
unit sphere thus gives a spherical harmonic of degree n, Y^(^). ■ 

Since the set {^n,i(C)}^i' n ^ i n (4-10) is a maximal linearly independent set 
of spherical harmonics of degree n, it serves as a basis for all such functions. 
We have just shown that Y n j(R£) is a spherical harmonic of degree n provided 
Y n ,•(£) is as well. Thus, we can write Y n j{R£) in terms of the basis functions: 

N(p,n) 

Y nJ (Rt)= >T C ej Y n/ (0- 
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Using this expression to rewrite the integral in (4.11), we can show that the 
matrix C defined in the above expression is orthogonal. Indeed, 



/N(p,n) \ /N(p,n) \ 

AT(p,n) 
. /i 1 ^ 



N(p,n) N(p,n) 

= E CkiCij5ke = E C\ k Ckj- 

k,£=l k=l 

For the following discussion, let £, 77 be two unit vectors. Let us consider 
the function given by 

N(p,n) 

Fn(tv)= E Y nM)YnM- ( 4 - 12 ) 
i=i 

Lemma 4.8 The function F n defined above is invariant under a rotation of 
coordinates. 

Proof Let R be a rotation matrix. Then, using the orthogonal matrix C 
discussed above, 

N(p,n) 

F n (Rt,R V )= E Y n,j(R^)Yn,j(Rv) 
0=1 

N(p,n) /N(p,n) \ / N(p,n) 

= E E C ^ Y nM) E CmjYnAn) 

j=l \ e=i J \ m=l 

so 

N(p,n) / N(p,n) 

F n (RZ,R V )= E ^(0^n, m (r?) E 

iV(p,n) 

= E ^(O^fa) = Fn&V), 

as claimed. ■ 

Since the dot product (£, ??) is also invariant under the rotation R, this 
suggests that F n (£, ?y) could be a function of (C,v) alone. In fact, this is the 
case , as shown in the following lemma. 
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Lemma 4.9 If F n is defined as in (4.12), then F n (£,rj) = p({£,T]}) where p(t) 
is a polynomial. 

Proof A rotation of coordinates leaves (£, rj) invariant and, by the above 
lemma, does not change F n (^,rj) either. For some —1 < t < 1, there exists a 
rotation R that sends 

^Ae' = (t,yr^,o,...,o), 

V ^rj' = (1,0,. ..,0). 

To see this, just rotate coordinates so that 77 points along the xi-axis. Then 
rotate coordinates around the xi-axis until the component of £ orthogonal to 
rj points along the X2-axis. 

Notice that (£,77) = (£',?/) = Since F n (£,rj) is a sum of products of 
spherical harmonics, which are each polynomials in the components of £ or 77, 
F n is a polynomial, say p, in the components of its arguments, i.e., 

F n (Z,r t ) = F n (e,r l ')=p(t,VT^). (4.13) 

We can impose another rotation of coordinates, again without changing F n or 
(£,77). Let R be the transformation that rotates vectors by ir radians about 
the xi-axis, i.e., 

e = (t, -7T^,o,..., 0), 

77' A 77" = (1,0,..., 0). 

Just as in (4.13), we conclude that 

Fn(ti, 7?) = F n (t", 77") = p(t, -VT^). (4.14) 

From (4.13) and (4.14), we see that yl — t 2 must appear only with even 
powers in F n . Thus, F n is really a polynomial in t and 1 — t 2 , which is just a 
polynomial in t. Since t = (£,77), the lemma is proved. ■ 

4.3 Legendre Polynomials 

Theorem 4.10 Let 77 = (1,0, . . . , 0) and let L n (x) be a harmonic homoge- 
neous polynomial of degree n satisfying 

(i) L n (rj) = 1, 

(ii) L n (Rx) = L n (x) for all rotation matrices R such that Rn = 77. 
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Then L n (x) is the only harmonic homogeneous polynomial of degree n obeying 
these properties. In particular, these two properties uniquely determine the 
corresponding spherical harmonic L n \ SP -i. Moreover, this spherical harmonic 
Ln(0 is a polynomial in (£,17). 

Proof Let £ G S p . Then there exist v £ S v ~ x such that (y,r\) = and 
t € [—1, 1] such that £ = trj + yl — t 2 v. Notice (£, rj) =t and that 

tx=t, $ + ---+e P = i-t 2 . (4.15) 

As in the proof of Theorem 4.4, we can write 

n 

L n {x) = y^x{hn-j(x2,X3, . . .,x p ), (4.16) 

J=l 

where the 7i n -j are homogeneous polynomials of degree n — j. Let R be a 
rotation matrix as described in property (ii) in the statement of the theorem. 
When R acts on x, it sends 



R i i 

X\ , X2 j • • • , Xp I r X\ , X2, • • • > Xp. 



By property (ii), 



— L n (x) — L n (Rx) — ^ ^ x j [/i n _jf(x2, • • • , Xp) — h n —j(x2, . . . ,x p 



3=1 



and using the linear independence of the x\, we see that all the h n -j(x2, ■ ■ ■ , x p ) 
are invariant under the rotation of coordinates R. Thus, these polynomials 

must depend only on the radius < lx\ H + ac|, i.e., 



/in-j(«2, •••,a;p) = c n _j (wa;|H h x£ , (4.17) 



where the c n -j are constants, and c n _j = for all odd n — j since the /i n _j 
must be polynomials. Property (i) gives us one of these coefficients, 

1 = L n (r]) = c 

since all the X2, . . . ,x p are zero for this vector. Since the form of the h n — j 
is given by (4.17), knowing cq is enough to determine the rest of the /i n _j 
using the recursive relation (4.6) and the fact that the c n -j are zero for odd 
n — j. Therefore, L n as well as the spherical harmonic L n \ SP -i are uniquely 
determined by properties (i) and (ii). Finally, using (4.16), (4.17), and (4.15), 

n 

j=0 
n—j=even 



60 



4. Spherical Harmonics in p Dimensions 



or 

n 

L n m))= E # (4.18) 

3=0 
n—j=even 

That is, L n (£) is a polynomial in t = (£,77), and the theorem is proved. ■ 

Definition We call the polynomial L n (C) introduced in Theorem 4.10 the 
Legendre polynomial? of degree n. Written in terms of the variable t, we 
denote it by P n {t). 

Remark We can quickly obtain some properties of Legendre polynomials. 
First, from (4.18), we can see P n (t) is a polynomial of degree n. We can also 
compute 

l = L n ( V ) = P n ((r,,r ] )) = P n (l). 

We can determine the parity of the Legendre polynomials in more than one 
way. First, from (4.18), we can see that if n is even, P n contains only even 
powers of t, while if n is odd, P n contains only odd powers of t. Alternatively, 
we can see that 



Pn(-t) = Pn((-C,V)) = L n (-0 = (-l)^n(O 

= (-i)»p n ((e,i/» = (-i)»p n (t). (4.i9) 

Either way, we determine that P n is even whenever n is even, and P n is odd 
whenever n is odd. 



In the following theorem, we demonstrate how to write the Legendre poly- 
nomials in terms of an orthonormal set of spherical harmonics. 



Theorem 4.11 (Addition Theorem for Legendre Polynomials) Let 

{^rijiO}^! ^ be an orthonormal set ofn-th degree spherical harmonics. Then 
the Legendre polynomial of degree n may be written as 

Pn{{M) = NM £ Y n,M)Yn,M- (4-20) 



2 We follow the naming convention used in [5] and common in physics. Mathematicians 
usually refer to these as Legendre polynomials only when p = 3 and as ultraspherical poly- 
nomials for arbitrary p. 
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Proof Since (4.20) is invariant under coordinate rotations, we can choose 
7] = (1, 0, ... , 0). Consider again the function F n (£, rj) defined in (4.12). Since 
we have already defined the vector rj, we will think of it as a fixed parameter in 
the function F n . For now, we will write F n {^\ rj) and think of F n as a function 
of one vector £. 

First, notice that F n , being a linear combination of spherical harmonics 
Yn,j(,0 (since we consider the Y n j(rj) to be constants), is itself a spherical 
harmonic, i.e., the restriction of a harmonic homogeneous polynomial to the 
unit sphere. Next, notice that any rotation of coordinates R that leaves rj 
fixed leaves the function F n invariant. Indeed, using Lemma 4.8 or Lemma 
4.9, 

F(R£; rj) = F(R£; Rrj) = F(£; rj) for all R such that Rrj = rj. 

Notice that in the application of the rotation R, we only rotate the coordinates 
of £ in F n since we consider rj a fixed parameter, though this made no difference 
in the calculation. 

Let us normalize the function F n by dividing it by the constant F n (rj;rj). 
Thus, we have found a function F n (£;rj)/F n (rj;rj) that obeys all the proper- 
ties of the Legendre polynomial described in the statement of Theorem 4.10. 
Since, by the same theorem, these properties uniquely define the Legendre 
polynomial, we conclude 

To complete the proof, we will compute F n (rj,rj). Since, by Lemma 4.9, 
F n (rj,rj) depends only on the inner product (rj,rj) = 1, it is a constant. Thus, 

F n (rj, 7/)n p _i = J F n (i],rj)d£l p - 1 

N(p,n) 



J E Y nt:j {rjf <m p ^ 

riesp- 1 • J=1 
£ J Y nJ (ri) 2 dn^ 1 =N(p,n), 



r/eSP- 1 

using the orthonormality of the n-th degree spherical harmonics. We see that 

p ft \ o N (p> n ) 



and we are finished. 
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In addition, we can expand spherical harmonics in terms of Legendre poly- 
nomials. To show this, we will need first the following result. 

Lemma 4.12 For any set {Y n j}j =1 of k < N(p,n) linearly independent n-th 
degree spherical harmonics, there exists a set {t]i}f =1 of unit vectors such that 
the k x k determinant 

Y n ,l(m) Yn,l(m) ■■■ Y n,l(Vk) 
Yn,2(m) Y n,2(m) ■■■ Y n,2(Vk) 

Y n ,k(m) Y n,k{m) ••■ Y n,k(m) 

is nonzero. 



(4.21) 



Proof We will prove the lemma by induction on k. First, consider a linearly 
independent set {^,1} of one spherical harmonic of degree n. Y n \ cannot be 
the zero function. Then, there exists a unit vector, call it 771, such that the 
determinant | ^,1(771)! = Y n \{j]\) 7^ 0. Thus, the lemma holds for the case 
k = 1. Now, suppose the lemma holds for some k = £ — 1 < N(p, n) — 1, and 
let {Y n j}j =1 be a set of I < N(p, n) linearly independent n-th degree spherical 

harmonics. By the induction hypothesis, we can choose a set }f=i of unit 
vectors such that the (I — 1) x {£ — 1) determinant 



A, 



Y n,i(m) 

^n,2(??ly 



Y nM 
^,2M 



Yn,2(V£-h 



7^0. 



(4.22) 



Now consider the spherical harmonic defined by the £ x £ determinant, 



A 



Y n,i(m) 

Y n,2(Vl) 

Y n,e-i(m) 
Y nAm) 



Y n.i(m) 

Y nM 

Y n,e-i(m) 

Y n,t{m) 



Y nt i(r)£-i) 

Y n,2 (W-l) 



Y n,l(0 
Y n,2(0 



Y n ^-i{r}i-\) Y n> i-i(£) 
Y nAm-i) Y nA0 



(4.23) 



If we compute this determinant by performing a cofactor expansion down the 
last column and indicating by Aj the minor determinant corresponding to 

e 

A = ^(-i)^A,y nj (0, 
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we can see that we will have a linear combination of spherical harmonics 
Y n ,j(0 which are linearly independent. Thus, for the determinant to vanish 
identically, all the coefficients of the Y n j(£) must vanish. But notice that the 
coefficient of Y n ^{£) is the determinant (4.22) which does not vanish. Thus, the 
spherical harmonic given by the determinant (4.23) is not the zero function; 
i.e., there exists a unit vector £ = rji such that the determinant (4.23) is 
nonzero. Therefore, the lemma holds for the case k = i and, by induction, for 
all k < N(p, n). U 

Theorem 4.13 For any spherical harmonic Y n (£) of degree n, there exist 
coefficients a k and unit vectors rj k such that 

N(p,n) 

Yn(0= £ a k P n ((£, Vk )). 

k=l 



Proof Let be an orthonormal set of n-th degree spherical har- 

monics, and let Y n (£) be any spherical harmonic of degree n. For a unit vector 
rj, we can write 

N(p,n) 



p-i 



N(p, n) 



E Y nA0Yn,M, 
3=1 



(4.24) 



by Theorem 4.11. Let us choose a set of unit vectors rjj such that the deter- 
minant (4.21) with k = N(p,n) is nonzero; this is possible by Lemma 4.12. 
Replacing rj by rjj in (4.24) for 1 < j < N(p, n) creates the system of equations 

/ ^«£,»7i» \ / Y nil { m ) ••• Y n)N(p>n) ( m ) \f Y nA (0 \ 

Y n ,i(m) ■■■ Y n 

,N(p,n) 



N(p, n) 



p-i 



Pn((Z,m)) 

\Pn((t,riN(p,n)))J 



Yn,2(0 
\Y n ,N(p,n) 



Vn,l('7iV(p ) n)) ••• Y n 

,N(p,n) (VN(p,n))J 

The determinant of the N(p, n) x N(p, n) coefficient matrix on the right-hand 
side of this system has nonzero determinant by our choice of the vectors rjj. 
Thus, this system is invertible; i.e., there exist coefficients C£ m such that 

N(p,ri) 

YnAO = £ c^ m P n (((, Vm )), for each 1 < I < N(p,n). (4.25) 

m=l 

Now since any spherical harmonic, and in particular Y n (t;), can be expanded 
in terms of the basis functions {Y ri: i(^)}^}^' n \ there exist coefficients a k such 
that 

N(p,n) 

k=i 
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by (4.25), and the theorem is proved. 



We will now list and prove several basic properties of spherical harmonics 
and Legendre polynomials, some of which are useful for making estimates. 
The proofs are straightforward. 



Lemma 4.14 For any spherical harmonic Y n (£), 



n 



p—i 



(4.26) 



Proof Let Y n (£) be any n-th degree spherical harmonic, and let {Y n j}^j^' n ^ 
be an orthonormal set of such functions. Then, we can expand Y n (rj) in terms 
of this basis; i.e., for some coefficients a,-, 

N(p,n) 

Yn{v) = a 3 Y n,j(.V)- 
3=1 

Using this expansion and Theorem 4.11 to rewrite Pn((£,,i]))y the right-hand 
side of (4.26) becomes 



N(p, n) 



N(p,n) 

a 3 Y n,j(v) 

3=1 



O p _i 
N{p, n) 



N(p,n) 



k=l 



dUp-i, 



or 



N(p,n) 

Y a j Y n,k(0 
j,k=l 



Y n> j(rj)Y ntk (ri) dQ } 



p-i 



N(p,n) 

]T a j Y n>j (0 = Y n (0, 
j'=i 



as required. 

Proposition 4.15 The Legendre polynomials P n (t) are bounded 

\Pn(t)\ < 1, for allte [0,1] 
and obey the following normalization condition 

tesp- 1 



N(p, n) 



(4.27) 



(4.28) 
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Proof We rewrite P n (t) 2 using Theorem 4.11. Then, we use the Cauchy- 
Schwarz inequality, viewing the sum of products in this expression as a dot 
product, to derive the required result: 



N(p,n) 



< 



n 



p—i 



N(p,n) 



N(p,n) 



a, 



N(p,n) 



so 



Pn«£, V)) 2 < Pn{(£, 0)Pn((v, v}) = ^n(l) 2 = 1, 

proving the first result. 

Again, we use Theorem 4.11 to rewrite the integral on the left side of (4.28) 

as 



/ 



which then becomes 



2 W(p,n) 

N(p, n) 2 



" 1 Y, Yn >i (O y n,j (v)Y n ,k (QYn,k (v) dSlp-i 
j,k=l 



q2 N(p,n) 



N(p, n 



^ P n ((v,V)) ~ ^ 



iV(p, n) 



iV(p, n) 



thus proving the second result. ■ 

Proposition 4.16 The spherical harmonics Y n (£) satisfy the following in- 
equality 



\Y„(0\ < 



N(p, n) 



Y n {r,) 2 dQp-L 



(4.29) 



rieSP- 1 

Proof We start by taking the square of equation (4.26): 



Yn(H) 



2 N{n,p) 2 



fl 2 , 
p—i 



Yn^Pnd^^dQp-! 
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Viewing the integral as an inner product, we apply the Cauchy-Schwarz in- 
equality, 



Yn(0 2 < 



N{n,pf 

n 2 ; 

p— 1 



y n (r/) 2 dfi p _i 



P„((C,7?)) 2 CKV-I 



Thus, 



S*>- 



where we used property (4.28) in the last step. ■ 

We will now begin to investigate the properties of Legendre polynomials as 
orthogonal polynomials. Let us rewrite the integral in (4.28). Note that the 
integrand depends only on the inner product (£, rj) and that we integrate £ over 
the surface of the (p— l)-sphere. We can take advantage of these observations 
to reduce the (p — l)-dimensional integral to a one-dimensional integral. 

Since we integrate over the entire sphere, we can perform any rotation 
of coordinates without changing the value of the integral. Let us impose a 
coordinate change R that aligns the unit vector r\ along the which we 

will picture pointing "north," i.e., take r\ = (0, . . . ,0, 1). If we let t = (£,77), 
we can write the unit vector £ as 



£ = {£, v)v + (£ - v)v) = tv + Vi-t 2 

for some unit vector v normal to r/. Notice that any such v gives the same 
value for the inner product (£, r\) and thus the same value for the integrand in 
(4.28). Note further that the collection of all such vectors v, 



\u\ = l,{u,r J ) = 0}, 



forms a parallel of the (p — l)-sphere which is a (p — 2)-sphere: 

{i/Gf: \v\ = l,v p = 0} = S P ~ 2 . 

A visual representation is depicted in Figure 4.1. 

To get some intuition, let's think of the familiar 3-dimensional case. The 
2-dimensional sphere can be thought as a sum of infinitesimal rings oriented 
along the parallels of the sphere. The infinitesimal ring defined by the az- 
imuthal angle 9 has an infinitesimal thickness d6 and, therefore, it corresponds 
to a solid angle 

dQ 2 = (27rsin6>) dO. 
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Figure 4.1: In the reduction of a spherically symmetric (p — 1) -dimensional integral to a 
1-dimensional integral, we imagine the (p — l)-dimensional sphere as a sum of infinitesimal 
(p — 2)-dimensional spheres (parallels) of thickness d8. 



The term inside the parenthesis is the length of the ring; multiplied by its 
thickness it gives the "area" of the ring. We can rewrite dQ.2 as —2ird(cos9) 
and thus 

dQ,2 = dt 

in terms of the variable t. 
Similarly, in p dimensions 

dn p _i = %- 2 (sm ey- 2 d6 . 

The factor (sin 8) p ~ 2 is easy to explain: The radius of the (p — 2)-sphere is 
R = sin (9; its volume will be proportional to RP~ 2 . It is straightforward to 
express dVt p -\ in terms of t: 

dQ p ^ = -n p - 2 (sm6) p ' 3 d(cos9) 

Now we are ready to write (4.28) as a 1-dimensional integral over t. Indeed, 

l 

J Pni^Ti)) 2 d%^ = j P n (t) 2 (l-t 2 ) P -^n p _ 2 dt, 
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which implies 

1 

*2/ n ,2^ A . _ ^p-i 



N(p,n)Q p - 2 ' 



This expression gives us the norm of the n-th Legendre polynomial with respect 



p-3 



to the weight w(t) = (1 — t ) 2 ; namely 



|P„(t)IU = V(Pn(t),Pn(t)) w = l/ ^r/^o ■ (4-30) 

Note that in coming up with this fact, we have essentially proved the 
following lemma, which will be used to show the next few results. 

Lemma 4.17 Let n be a unit vector and f be a function. Then 

1 

f((tv))dn p ^ = ft p _2 J f(W - t 2 ) {p -' i)/2 dt. 

tesp- 1 -1 



Theorem 4.18 Any two distinct Legendre polynomials P n (t),P m (t) are or- 
thogonal over the interval [—1, 1] with respect to the weight (1 — t )~^~ ■ That 
is, 

1 

J P n (t)P m (t)(l - t 2 )^ dt = 0, for n^m. 
-l 



Proof Let n = (1,0, .. . ,0). The Legendre polynomial P n ({£,v)) ^ s equal to 
the spherical harmonic L n {^) having the properties listed in Theorem 4.10. 
By Theorem 4.6, 

0= J L n {i)L m {i)d%^= j P n ((e,r?))P m ((e,r?))^ P -i, for n ± m. 

Since the integral on the right-hand side of the above equation depends only 
on the inner product (£, 77), we can use Lemma 4.17 to rewrite this integral as 

1 

fip-2 J P n (t)P m (t)(l - t 2 ) 2 ^ dt = 0, for n^m, 
-1 

thus completing the proof. ■ 
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The above theorem allows us to use the results from Chapter 3 to write 
down more properties of the Legendre polynomials. Using the comment just 
above Section 3.4, we can find the Rodrigues formula for the Legendre poly- 
nomials in p dimensions. Theorem 4.18 tells us that the Legendre polynomials 
are orthogonal with respect to w(t) = (1 — i 2 )( p_3 )/ 2 so that 



P n (t) = c n {l-t 2 )-^' 2 



Cn(l - t 



2x(3-p)/2 



d_ 
dt 

d_\ n 
dt 



:i _ t 2 ) ( P -3)/2 (1 _ t 2 ) n 



l-t 



2\ra+(p-3)/2 



and it remains to compute c n . We know that -P n (l) = L so 

1 = Cn (l _ t 2 )( 3- P )/2 fd\ n (1 _ t 2 r+ ( P -3)/2 



Carrying out two differentiations, 



l=cjl-t 



2n(3- P )/2 



n-1 



On.il ~ t 



2n(3- P )/2 



dt 

dt J 



n-2 



n + (p - 3) 



n + (p - 3) 



-2t)(l -t 



t=i 



2Nn-l+(p-3)/2 



(-2t) 2 (l - t 



2\n-2+(p-3)/2 



+ ••■ 



t=l 



where we have used the falling factorial notation and left terms that are higher 
order in 1 — t 2 in the " • • • " because they will vanish when we substitute t = 1. 
Continuing the pattern, 



1 = Cn (l - t 2 )( 3 -P)/ 2 Un + (p- 3)/2) n (-2t) n (l - t 2 )"-«+(p-3)/ 2 + . . 
= c n (n+(p-3)/2) n (-2)", 

implying that 



t=i 



(-ir 



2 n („+(p_ 3 )/2) n - 



We have shown the following. 



Proposition 4.19 (Rodrigues Formula for Legendre Polynomials) 

^n+(p-3)/2^ (4.31) 



p ( t ) - ( 1 )" f1 _ ^(3-p)/2 M V (-r _ j2\t 

^-2«(n + (p-3)/2) n (1 [dt) (1 ^ 



The above Rodrigues formula allows us to prove the following two proper- 
ties of the Legendre polynomials. 
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Proposition 4.20 In p dimensions, the Legendre polynomial P n (t) of degree 
n satisfies the differential equation 

(1 - t 2 )P^(t) + (1 - p)tP' n {t) +n(n+p- 2)P n {t) = 0. (4.32) 

In what follows, we will often write P n and w instead of P n {t) and w(t) for 
simplicity. A portion of the proof will be left as a straightforward computation 
for the reader. 

Proof Consider the expression ^ can use the product 

rule to rewrite this as 

(1 - t 2 )^l 2 {-2t)P' n + (1 - t 2 )^' 2 Pl 

or, 

w [{l-p)tP' n + {l-t 2 )P^\. (4.33) 

The term in brackets is a polynomial of degree n which can be written as a 
linear combination of the first n Legendre polynomials. Hence 

j n 

U{l-t 2 )wP' n ] =wJ2c j P j - (4.34) 

3=0 

Multiplying each side by Pk, where < k < n and integrating over the interval 
[-1, 1], we get 

l 

c fc ||P fc || 2 = J P k j t [{l-t 2 )wP' n ] dt. 

-l 

Integrating the right-hand side by parts now twice and noticing that the 
boundary terms vanish, we find 

i 

c fe ||P fe || 2 = P fc (l-t 2 )^ti- / P' k (l - t 2 )wP' n dt 
* ' i 

_i 

1 

= -P^(l-t>Pn|^ + j {j t [P k {l-t 2 )w]\p n dt. 

The expression in curly brackets can be written as wq k , where q k is a poly- 
nomial of degree k. We can then use Proposition 3.3 for the two integrals to 
discover that c k = for all k < n. Returning to equation (4.34), this result 
implies that 

|[(l-t 2 X] =wc n P n . 
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To determine c n , we will compute the coefficient of the highest power of t 
on each side of the above equation, using Newton's binomial expansion. By 
keeping only the highest-order term in each binomial expansion, the reader 
can show that the left-hand side becomes 

2 n {n+{p _ 3)/2)n ^ n + P W + > 

while the right-hand side becomes 

(_1)(p-3)/2 {2n + p-3) n +n+p _. 



2 n (n + (p- 3)/2) r 



+ 



Comparing the above equations, we conclude that c n = — n(n+p— 2). Inserting 
this into (4.34) and using (4.33), we find 

(1 - t 2 )P' n ' + (1 " V)tP'n + n(n + p- 2)P n = 0, 

completing the proof. ■ 

Proposition 4.21 The Legendre polynomials in p dimensions satisfy the re- 
currence relation 

{n + p- 2)P n+l {t) - (2n + p - 2)tP n {t) + nP„_i (t) = 0. 

Just as in the previous proof, we will leave part of the following proof to 
the reader. 



Proof We know from Proposition 3.5 that a relation of the form 

P n+l - {A n t + B n )P n + C n P n ~x = (4.35) 

exists. We can quickly determine B n . Recall from (4.19) that P n is even 
(respectively, odd) whenever n is even (respectively, odd) . Rewriting the above 
equation as 

Pn+l A n tP n + C n P n —i — B n P ni 

we have an odd polynomial equal to an even polynomial, unless both sides 
vanish. Since the first option is not possible, we conclude that both sides of 
the equation must cancel, which implies that B n = 0. We will use the notation 
and results of Proposition 3.5 to compute A n ,C n . Keeping only the highest- 
order terms in each binomial expansion in the Rodrigues formula (4.31) and 
carrying out the derivatives, it is straightforward to show that the leading 
coefficient of P n is 

(2n+p-3) n 
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from which we determine the leading coefficient of P n +i-, 

(2n + p- l) n +l 



K+1 2"+i(n + (p-l)/2) n+1 



allowing us to compute 



Now, using (4.30), 



_ fc n+ i _ 2n + p - 2 
k n n + p- 2 



^ An 1 1 P<, 

'in-1 ll-fn-1 



2 



n || - 1 n || 

4 7 Tip 77|2 



2n+p — 2 n+p — 3 iV(p, n — l)f2 p _2 

n + p-2 2n + p-4 iV(p,n)f]p_2 

Using Theorem 4.4, we can compute 

c n 



n + p — 2 

Inserting these results into (4.35) and multiplying by n + p — 2, we find the 
required result. ■ 

It is also interesting to observe that once we find the Legendre polynomials 
in p = 2 and p = 3 dimensions, they are determined for all higher dimensions. 
In particular, the following theorem proves, more specifically, that the Legen- 
dre polynomials for even p can be found from the p = 2 case and that those 
for odd p can be found from the p = 3 case. For the next theorem, we will let 
Pn,p(t) denote the n-th degree Legendre polynomial in p dimensions. 

Theorem 4.22 For all j = 0, 1, . . . , n, 

((p-l)/2)j (aiy 



Pn-J,2 j+P (t) - ( _ n) . (n+p _ 2) . [j t ) P nAt)- 
Proof Differentiating (4.32) once with respect to t, we get 

(1 - t 2 )P",' P + ("I - P)tP'n, P + {n-l)(n + p- l)P; p 
but this is just (4.32) again with the following substitutions: 

p i — ► p + 2, 
n i — > n — 1. 
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Thus, solving this new differential equation for P' will give 

K^tt) « Pn-l,p+2(t). 

Continuing in this way, if we differentiate (4.32) j times we see that 

^ 3 



j£ ; Pn,p{t) OC P n -j iP+ 2j(t) 

for any < j < n. Since i- > n _j iP+ 2j(l) = 1, in order to create an equality we 
must divide the left side of the above equation by its value at t = 1, i.e., 



Pn-j,p+2j (t) 



d ^Pn,p(t) 
d \3 



* dt , 



Using the Rodrigues formula (4.31), 

d \ j o /.M _(-n) j (n + p-2) j 



dt) Pn ^\t=i- ((p-l)/2)j 
as the reader can verify by computation. Therefore, 

((p-l)/2), fd\ 



Pn-j,2j+p{t) — ~, r~ ; ; -r- -77 Pn,p(t) , 

{-n) j (n + p-2) j \dtj 
completing the proof. ■ 

The following lemma will be useful in the proof of the theorem that follows 
immediately afterwards. 

Lemma 4.23 Let £ and £ be unit vectors, f a function, and -F((,,£) given by 



Then, 

1 

F(c,e) = ^- 2 Pn((e,o) J f(t)p n (t)(i-t 2 )^/ 2 dt. 



Proof First, notice that F is invariant under any coordinate rotation R, i.e., 
F{R(,R£) = F(C,0- Indeed, the rotation 
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can be undone by a change of variables 

r\ i — \ rj = Rrj, 

in the integration. Thus, we are allowed to choose our coordinates such that 

t = (1,0,. ..,0), c = (^,v / i-s 2 ,o,...,o), 

so that (C, £) = s. 

Let us think of F as a function of £ alone and £ as a fixed parameter; we 
will write .F(£;£). The argument of F, namely £, only shows up inside the 
Legendre polynomial P n ({n, 0), which is a spherical harmonic — in particular, 
a harmonic homogeneous polynomial of degree n in the variables £i, £2, • • ■ , Cp- 
Therefore, the function F is also a harmonic homogeneous polynomial in the 
components of £, i.e., in s, Vl - s 2 . But, since we equally well could have 
chosen coordinates such that Q = (s, — v 1 — s 2 , 0, . . . , 0), F must really be a 
polynomial in s. 

Then, being only a function of the inner product s = (C>£)> we see that 
F satisfies all the defining properties of the Legendre polynomials in Theorem 
4.10 except (i). We conclude 

_F(C,£) = cP n (s), for some constant c. 

We can determine the constant by considering the case s = 1, i.e., £ = £, 
where -F(£,£) = cP n (l) = c. Using Lemma 4.17, 

1 

C = J f((c,v))Pn({v,0)dn p ^ = n^ 2 j f(t)p n (t)(i-t 2 )^/ 2 dt, 

as sought. ■ 

Theorem 4.24 (Hecke-Funk Theorem) Let £ be a unit vector, f a function, 
and Y n a spherical harmonic. Then, 

1 

/ /((£, v))Y n (v) dOp-l = np_ 2 y„(0 / - *) (P_3)/2 dt. (4.36) 

Proof By Theorem 4.13, there exist a set of coefficients {afc}^ P ' n ^ and a set 
of unit vectors {Cfc}fe=i such that 

N(p,n) 

Y n{i)= ^2 a k P n((V,Ck})- 
k=l 
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Then, we can rewrite the left-hand side of (4.36) as 



N(p,n) 




k=l 



or, using Lemma 4.23 




N(p,n) 



£ a k Pn{{Uk))\ f(t)P n (t)(l-t 2 ) ip - 3)/2 dt 




That is 



i 




-i 



which is exactly the right-hand side of (4.36). 



We wish to find now an integral representation of the Legendre polynomi- 
als. Towards this goal, we first prove a lemma. 

Consider a vector n of the (p — l)-dimensional sphere S^ 1 . Without loss 
of generality, we may take it along the xi-axis. With {n} 1 - we indicate the 
set of all vectors that are perpendicular to rj; this is obviously a hyperplane. 
The set {n} 1 - n S^ 1 is then the equator of S^ 1 (a (p — 2)-sphere) that is 
orthogonal to n; we will indicate it by S%~ . 

Lemma 4.25 Let rj = (1, 0, . . . , 0) and x £ W . Then the function 



when restricted to the sphere, is the n-th Legendre polynomial: L n {^) = P n ({£,, ?]))■ 

Proof Clearly, L n (x) is a polynomial in the components of x. By the bino- 
mial theorem, this polynomial is homogeneous of degree n. We can also see 
that L n {x) is harmonic. Indeed, applying the Laplace operator on (4.37) and 
switching the order of differentiation and integration, the integrand becomes 




(4.37) 
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or 

p 



or 



But 



d r 



l>i + ^') 2 = fa> ^ " (C. + 2i(77, o = 0, 



since {t?,C} is an orthonormal set. Also, notice 



dJ7p_ 2 = 1. 



Now, let i? be any coordinate rotation leaving rj invariant, i.e., let R be any 
rotation about the xi-axis. The integrand of L n {Rx) is then 

[(Rx, V ) + i(Rx, Of = [<z, »/) + *<z, #0] n , 
which can be reset to 

[(x,rj) + i(x,()] n 

by a change of variables ( — > £' = RQ and thus L n (x) is invariant under all 
such coordinate rotations. 

We have shown that L n {x) has all the properties described in Theorem 
4.10, so that, when restricted to the sphere, it becomes the re-th Legendre 
polynomial. ■ 

Now we can prove the following integral representation for P n (t). 

Theorem 4.26 

l 



n 



p-2 



-1 
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Proof From Lemma 4.37, we have 



Choose a constant t and unit vector v such that £ = tn + \/ 1 — t 2 v and (i/, 77) = 
0. Then, the above equation becomes, using Lemma 4.17 and replacing p by 
p-1, 

Pn(t) = fT~ 2 j + -* 2 <^,C>)" d%-2 
l 

ft 



o 

as sought. ■ 

4.4 Boundary Value Problems 

We conclude this discussion with an application of the ideas we have devel- 
oped to boundary value problems, where they display most of their physical 
importance. 

We know from Proposition 3.15 that if a set of functions in a Hilbert 
space is closed, then it is complete. In the following theorem, we will see 
that a maximal linearly independent set of spherical harmonics of all degrees 
is closed and thus complete. This result allows us to develop expansions of 
functions as linear combinations of spherical harmonics, which will be useful 
in application to boundary value problems. In what follows, let 

S = {Y n j : n G No, 1 < j < N(p, n)}, 

be a maximal set of orthogonal spherical harmonics, and let ^ n • denote the 

Theorem 4.27 Let the function f : 5 P_1 — > R be continuous. If f is orthog- 
onal to the set S, i.e., if 

J f(0Ynj(0 dfip-i = 0, for all n,j 

then f is the zero function, i.e., 

/(C) = 0, for allies?- 1 . 
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The requirement that / be continuous is actually too strict, and the theo- 
rem really applies to all square-integrable functions /, i.e., all / such that 



J /(0 2 <ffip-i<°o. 



However, we will only prove the weaker version of this theorem. 

Proof We will prove it by contradiction. Suppose that / satisfies the hy- 
potheses of the above theorem, i.e., is continuous and orthogonal to S, but is 
not the zero function. Then there exists some rj G S p such that f(r)) 7^ 0. 
We can assume that f(rj) > 0, for if f{rj) is negative we could consider — / 
instead. By the continuity of /, there is some neighborhood around rj on the 
sphere where / is positive. That is, there exists a constant s such that /(£) > 
whenever s < (£, r/) < 1. Define the function 

I if - 1 < t < s. 

For s < {£,rj) < 1, the product /(£)?/>((£, rj)) is positive, and it vanishes for all 
other £. Thus, 

j /(^((^»ffl H c > ' ( 4 - 38 ) 

For simplicity, we will indicate the above integral by c. By the Weierstrass 
approximation theorem (Proposition 3.10), we can find a polynomial p(t) for 
any given e > such that 

\ip(t) - p(t)\ < e, for all t G [-1,1]. 

For any such e and p(t), 



/(0 [^«e,»7»-p«e,»7»] ^ P -i< 



/(C) [V'((C^))-p((C,r 7 ))] dfip-i 
< J \f(£)\\il>((£,ri))-p((S,v))\dtop-i 



Since /(£) is continuous and £ G S* 1 '" 1 , there exists an M such that M > |/(£)| 
for any £ G S ,p_1 . This implies 



/(0K&»7>)<KV-i>c-AfeJVl- 
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And since we can choose e arbitrarily small, 

J /(£M<M)<KVi>0- ( 4 - 39 ) 

For the remainder of the proof, we fix e and the corresponding p(t) for which 
this expression is true. 

Now, let m denote the degree of p(t). We can write the polynomial p(t) as 
a linear combination of the first m Legendre polynomials since each P n (t) is 
of degree n; i.e., we can find such that 

m 

P (t) = Y J c k p k (t). 

k=0 

We can thus rewrite the integral in (4.39) as 

» m 

J /K)5]^((^»<«p-i. 

But since the Legendre polynomials are just a special collection of spherical 
harmonics in £, this integral must vanish by hypothesis, contradicting the 
assertion in (4.39). Therefore, our initial assumption that / is not the zero 
function must be false. ■ 



So for any reasonable function / defined on the sphere, we can write 

m)=J2^ Y n,M (4-40) 



To find c n >j>, multiply both sides of this equation by Y n 'y and integrate over 
the sphere. Using the orthonormality of the set S, we find 



C n ',j' 



J /(£)*V,i'(£)<ftVi- ( 4 - 41 ) 



5p-i 

This expansion will be used in the following demonstration. 



Problem Consider the following boundary-value problem. Find V in the 
closed unit ball B p , such that 



A p V = 0, and V = /(£), for all £ G 



(4.42) 
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Solution We know that harmonic homogeneous polynomials are solutions to 
the Laplace equation, as well as any linear combination of them. We have 
also seen in (4.7) that we can write each of these polynomials as a power of 
the radius multiplied by a spherical harmonic. Thus, we can construct the 
solution to (4.42) as a linear combination of r n Y n j((f) terms. To satisfy the 
boundary condition, we can use the coefficients in (4.41) to find 

^ = E rn W^) = E r " F ^) / f(v)YnAv)dn p ^, (4.43) 
n,j n,j 5p _ 1 

thus solving the problem with the solution being in the form of a series. ■ 

In the next subsection, we will learn to solve this problem by another 
method. Since the solution of this boundary- value problem is unique (as we 
know from the theory of differential equations), we can equate the answers. 
And, amazingly, this procedure will give us a generating function for the Leg- 
endre polynomials P n {t). 

Green's Functions 

Given a differential equation 

D x {f(x)) = , 

where D x is a differential operator acting on the unknown function f(x), the 
corresponding Green's function G is defined by the equation 

D x (G(x)) = 5(x - x ) , 

where the function 5 appearing in the right-hand side is the Dirac delta func- 
tion defined by 

5(x — xq) = 0, for all x / xq, 

and 

J 5(x - x ) d p x = 1, for all e > 0. 

Let's assume that the given differential operator is the Laplacian in p 
dimensions, that is, we seek the Green's function which satisfies 

A p G = 5(x-x ). (4.44) 

If we think of G as electric potential, then (4.44) describes the electric potential 
caused by a point charge at xq 
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Since the Laplacian is invariant under translations, we see that G can only 
depend on the distance from xo, i.e., on p = \x — xq\. When p / 0, G must 
satisfy Laplace's equation: 

= A p G(p) = ^G(p) 

i=l 1 



S^_d_( 9G(p) dp_ 
4^ dxi \ dp dxi 



i=i i=i 1 

= G"{ P ) + P —G\p). 
P 

We can easily solve this differential equation by separation of variables to find 

G(p) = ap + b, if p = 1, 

and 

if p > 2. This, in turn, implies 

G(p) = a In p, ii p = 2, 

and 

pP z 

We will focus on the last case but the reader should explore the cases p = 
1, 2 on his or her own to get a better understanding. To find the undetermined 
constant, we integrate the defining equation over a ball of radius e centered at 
the point xq- 



J A p Gd p x= J 8(x-x )d p x. 



Bf(x ) Bf(x ) 



By the properties of the Dirac delta function, the right-hand side is 1. The 
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left-hand side can be rewritten by the use of the divergence theorem: 
J A p Gd p x = J V P G •f(e p - 1 dn p -i) 

S? _1 («o) 
= a (2 -p) J2 p _i. 

Hence 



(2-p) 

Of course, differential equations come with boundary conditions. So, let's 
modify the previous problem as follows. Let's seek the Green's function which 
satisfies the same equation 

A p G = 5(x-x ), (4.45) 

for all x 6 B p (0) and subjected to the boundary condition 

G(0 = 0, for all £ G S*- 1 . (4.46) 

Equation (4.45) now describes the electric potential caused by a point charge 
at xq and an ideal conducting sphere with center at the origin. 
To construct G, we write 

G(x; Xo ) = G(p) + g = {2 _ p) ^_ lpP -2 + S, 

and require that g is harmonic in B p 

A p g = 

and cancel G on the boundary of B p : 

In fact, the functional expression of g is identical to that of G. However, there 
are two parameters that have to be fixed: the location of the singular point 
representing a point charge and the strength of the charge. The location of 
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the singular point of g cannot be inside B p (0). We will place the charge at the 
symmetric point x' to xo with respect to the sphere, i.e., the point x' that 
lies on the line passing through the origin and xq with \x' \ = l/|xo|- We see 
then that 

9{P>) " (2-p)^"- 2 ' 

where p' = \x — x' \. We must choose the charge to be of the correct strength 
so that G vanishes on the sphere. We easily see that the correct choice of g is 

g(p') = ' 



(2-p)f] p _ 1 (|x |p') p - 2 ' 
so that 



Indeed, (4.47) satisfies Laplace's equation. Using the law of cosines and letting 
be the angle between the ray from the origin to x and the ray from the origin 

to Xq, 



p =\j |x| 2 + |xq| 2 — 2|x| |xq| cos#, 



p =\l\x\ 2 + Ix'J 2 — 2\x\ \x'r,\ cosO = \ l\x\ 2 + | — — 2- — - cos 9, 

we see that (4.47) vanishes on the unit sphere, i.e., when |x| = 1. The method 
of constructing G as described above is known as the method of images. Physi- 
cists use it routinely without paying attention to the mathematical details! 
The reason is that if a solution is found for a boundary problem, it must be 
unique. 

Let's now return to the problem stated on page 79 and present an al- 
ternative solution using the results on the Green's function for the Laplace 
equation. 

Alternative Solution Green's theorem (as shown in the footnote of page 
55) for the function G and V, 

j (VA P G - GA p V) d p x = j (v^-G^j dQ^t, 

Bp SP" 1 

reduces to 



/ 



/BC 
V—dQp-!, 



Bp SP-i 
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since V is harmonic and G obeys (4.45) and (4.46). Since the left-hand side 
of the above equation becomes 

V8(x - xq) d p x = V(x ) J 5{x - x Q ) = V(x ), 

Bp Bp 

and since 



dG dG 



1 - \xq\ 



d£ d\x\ \ x \=i O p _i(l + |xo| 2 -2|2;o|cos6')P/ 2 ' 
we can write the solution, 

/ ^(TTSra^^- 1 - (4 ' 48) 

tesp- 1 

Obviously, this gives the potential as an integral representation. ■ 
Equating the two solutions, we can arrive at the following result. 

Theorem 4.28 In W we have 

00 i 2 
V r n N(p,n)P n (t) = — — — — — — jt . 

Proof Let xq = \xq\t]. Starting from equation (4.43), 



n,j teSP- 1 

oo . N(p,n) 

Ei^i" / E ^(o^,i^)^ P - 



n=0 tiesp- 1 1 i 



we rewrite it using Theorem 4.11, 

y(* ) = f> r / /(0^^iM&»7>)<aV-i 



n_0 £eSP- 



1 



fj p _i 



/ (e) E i x ° i n ) p " ( cos ) dJ v 1 • 



Since the function / is arbitrary, we can compare the above equation with 
(4.48) and set r = \xq\ and t = cosQ to complete the proof. ■ 
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This concludes our development of spherical harmonics in p dimensions. 
We have briefly considered an application to boundary value problems; we 
will not delve further into applications. We have achieved our main goal 
to study the theory of spherical harmonics and the corresponding Legendre 
polynomials in W. We urge the reader to seek out applications on his or her 
own. Perhaps search for instances in physics where spherical harmonics are 
used in M 3 and try to generalize to p dimensions. One could start with the 
multipole expansion of an electrostatic field (see [6], [7]) or the wave function 
of an electron in a hydrogenic atom (see [11], [8]). 



Bibliography 



[1] G. Arfken, H. Weber, and F. Harris. Mathematical Methods for Physicists, 
Fifth Edition. Academic Press, 2000. 

[2] F. Byron and R. Fuller. Mathematics of Classical and Quantum Physics. 
Dover Publications, 1992. 

[3] G. Folland. Fourier Analysis and its Applications. American Mathemat- 
ical Society, 2009. 

[4] S. Friedberg, A. Insel, and L. Spence. Linear Algebra. Prentice Hall, 
2002. 

[5] H. Hochstadt. The Functions of Mathematical Physics. Dover Publica- 
tions, 1987. 

[6] J. Jackson. Classical Electrodynamics, Third Edition. Wiley, 1998. 

[7] L. Landau and E. Lifshitz. The Classical Theory of Fields, Fourth Edition. 
Butterworth-Heinemann, 1980. 

[8] L. Landau and E. Lifshitz. Quantum Mechanics: Non-Relativistic Theory, 
Third Edition. Butterworth-Heinemann, 1981. 

[9] H. Royden and P. Fitzpatrick. Real Analysis. Prentice Hall, 2010. 

[10] W. Rudin. Principles of Mathematical Analysis. McGraw-Hill, 1976. 

[11] R. Shankar. Principles of Quantum Mechanics. Springer, 1994. 

[12] J. Stewart. Calculus: Early Transcendentals. Brooks Cole, 2007. 

[13] W. Strauss. Partial Differential Equations: An Introduction. Wiley, 2007. 



87 



